Technical Report RSC-37 


IHFBSBED SMSOES irElMf i P OBSp E " ;il ' CE ° P 
1-F.^chnson (Texas 


CSCL 14B 


G3/14 


W72-32472 


Unclas 

43478 


ON THE PERFORMANCE OF INFRARED SENSORS 
IN EARTH OBSERVATIONS 


by 


L. F. JOHNSON 



Reproduced by ' ~ — 

national technical 
INFORMATION SERVICE 

u S Deportment of Commerce 
Springfield VA 22151 



Technical Report RSC-37 


ON THE PERFORMANCE OF INFRARED SENSORS 
IN EARTH OBSERVATIONS 
by 

L. F. Johnson 


August 1972 


supported by 

National Aeronautics and Space Administration 
NASA Grant NsG 239-62 
NASA Contract NAS 9-11155 


( 



ABSTRACT 


The performance of infrared sensing systems is depen- 
dent upon the radiative properties of targets in addition 
to constraints imposed by system components. The unclas- 
sified state-of-the-art of infrared system performance 
figures is reviewed to indicate the relevance to system 
performance of target radiative properties. A theory 
of rough surface scattering is developed which allows the 
formulation of the reflective characteristics of extended 
targets. The thermal radiation emission from extended 
targets is formulated on the basis of internal radiation 
characteristics of natural materials and the transmissive 
scattering effects at the surface. Finally, the total 
radiative characteristics may be expressed as functions 
of material properties and incident and received directions, 
although the expressions are extremely complex functions 
and do not account for the effects of shadowing or mul- 
tiple scattering. It is believed that the theory may be 
extended to include these effects and to incorporate the 
local radii of curvature of the surface. 
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CHAPTER I 
INTRODUCTION 

The measurement or identification of some property of 
an object (target) without making physical contact between 
the target and the sensor is the essence of the remote 
sensing mission. Detection and measurement of any phenom- 
enon which exhibits a dependency on the target character- 
istics may serve to accomplish this mission, although elec- 
tromagnetic sensing is particularly suited to the remote 
sensing task. In fact, the predominance of electromagnetic 
remote sensing has become such that "remote sensing" effec- 
tively implies the detection and measurement of electromag- 
netic radiation. 

Remote Sensing 

The target characteristics which may influence elec- 
tromagnetic radiation (type of material, temperature, sur- 
face state) are geophysical in nature. As a result, the 
primary applications for remote sensing include identifica- 
tion or measurement of geophysical phenomena. Secondarily, 
the geophysical state may be related to various other 

/ 
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phenomena which are not directly detectable, but the basic 
problems involved in geophysical interpretation remain. 
Various non-geophysical situations which are accompanied 
by remotely detectable geophysical phenomena are conse- 
quently detectable by remote sensing techniques. For ex- 
ample, seed germination may be detected by measurement of 
the associated temperature change which is known to occur 
during germination. 

Remote sensing using airborne sensors has made possi- 
ble the rapid acquisition of geoscience data over large 
areas without significant increase in cost over conven- 
tional methods. The shorter acquisition time associated 
with airborne coverage is responsible for a consequent im- 
provement in accuracy of the data. Temporal variations 
are reduced to a minimum and the cost, although higher per 
unit time, is comparable due to the time reduction. 

Successful application of remote sensing techniques 
has been made in numerous fields of natural and physical 
science: meteorology, oceanography, geology, astronomy, 

agriculture, ..., our environment in general, offer a mul- 
titude of potential applications. Remote sensing has been 
used effectively in monitoring environmental pollution, 
thermal mapping of the Earth's surface, ground water pros- 
pecting, detection and location of forest fires and sub- 
surface coal fires, and early detection of diseased crops. 



The primary limitation on remote sensing occurs in 
data interpretation. Hardware systems are capable of 
accurate measurements to a degree that far surpasses human 
or automated analysis ability. Failure to consider cer- 
tain aspects of the remote sensing phenomenon and applica- 
tion of simplifying assumptions are responsible for ques- 
tionably valid conclusions about experimental data. 

An excellent introduction to remote sensing and its 
applications is given by D. C. Parker, et. al. [1], [2], 

and reports of recent research are available in the Pro- 
ceedings of the Uni vers ity of Michigan International Sym- 
pos iums on Remote Sens ing of Environment [ 3-9 ] . 

Radiation Phenomena 

The remote sensing capability afforded by electromag- 
netic radiation may be attributed to three recognizable 
phenomena: emission, reflection, and propagation. Every 

object having nonzero absolute temperature and emissivity, 
emits and reflects electromagnetic energy which exhibits 
a dependency on various characteristics of the target. 

The radiation travels (propagates) away from the target 
in all directions and becomes incident upon the remote 
sensor. The frequency, phase, direction, amplitude, and 
polarization of the radiation indicate such target proper- 
ties as temperature, type of material, and surface state 
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(roughness). Consequently, information about the target 
is conveyed to the sensor via the transmitted energy. 
Theoretically, the target characteristics may then be de- 
termined through proper measurement and interpretation of 
the resulting data, i.e. by remote sensing. 

The remote sensing situation, however, is complicated 
by several physical phenomena. If the particular property 
of interest produces a measurable influence on the radia- 
tion energy, remote sensing is possible; otherwise it is 
not. Clearly, the influence is measurable if it is not 
critically deteriorated by influences from other sources. 
The total radiation received by the sensor is composed of 
radiation originating from various sources and is conse- 
quently dependent upon the characteristics of other objects 
in addition to the target. Radiation from non-targets may 
reach the sensor directly or by reflection from the target 
itself. The total effect of this background radiation 
contribution is to obscure the target radiation character- 
istics. Also, several characteristics of the target may 
affect the radiation in a similar manner such that the 
radiation characteristics may not be unambiguously related 
to the target characteristics. The total radiation, 
emitted and reflected from the target is also perturbed as 
it passes through the intervening medium to reach the sen- 
sor; scattering and absorption add further complexity to 
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the problem. The resulting heterogeneous combination of 
energy contributions and radiation phenomena does contain 
information, although often obscure, about the target char- 
acteristics. Interpretation of the data in view of the 
total radiation situation is the difficult problem of the 
data interpreter. 

Data Interpretation 

( 

The remote sensing capability is hindered by the 
existing inability to interpret the resulting data accu- 
rately, frequently as a result of the interdisciplinary 
nature of remote sensing. Several factors contribute to 
the problem of interpretation, including the performance 
of the measuring instrument, the geophysical situation, 
and the particular application. 

The increased use of remote sensors is accompanied by 
an increased misuse due to a lack of knowledge about the 
instrument capabilities and limitations. The ability of a 
remote sensing instrument to identify characteristics of 
an object is necessarily dependent upon the properties of 
the components which comprise the measuring device. 

Failure to consider these "sensory" constraints, some of 
which are common to all systems, while others vary with 
information requirements and system design, seriously 
impairs the validity of experimental results. 



Much of the data interpretation problem stems from a 
lack of thorough understanding of the extra-sensory physi- 
cal phenomena involved. The exact nature of various phe- 
nomena important in electromagnetic sensing, including 
natural radiation emission, scattering of electromagnetic 
waves from rough surfaces, and atmospheric transmission of 
electromagnetic radiation are not understood in detail. 

The extra-sensory interpretation problem may be sim- 
plified by any of a number of assumptions about certain 
aspects of the total radiation situation. Any of the tar- 
get, background, or intervening medium effects may be 
neglected or assumed to be of a certain nature at the cost 
of some accuracy., A common assumption is that the radia- 
tion characteristics of the target are those of a black- 
body, although in nature objects may approach but do not 
possess blackbody characteristics. The state-of-the-art 
is such that these assumptions are necessary to obtain 
results that are valid within the restrictions imposed by 
the assumptions. 

The interpretation problem may also be simplified for 
certain applications in which the reflected component of 
the received radiation contains the desired information. 
The solution is to overwhelm all of the natural radiation 
components with a generated radiation component having 
specific characteristics. After reflection from the tar- 
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get, the radiation may be remotely sensed as before, but 
interpreted more easily. Radar, for example, utilizes 
this approach. The complex combination of natural radia- 
tion components is suppressed as noise and the single 
generated component is the signal to be detected and in- 
terpreted. 

Generation of this additional component is the dis- 
tinguishing characteristic of an active sensing system as 
compared to a passive system which is designed to detect 
naturally occurring radiation. The utility of active and 
passive systems is dependent upon the particular applica- 
tion, or more specifically, on the particular geophysical 
situation to be detected. For example, active systems are 
useful in measurements of surface roughness or target-to- 
detector distance as in radar, however, passive systems 
prove more valuable for measuring material properties. 

The nature of the applicability of each type of system is 
revealed by comparison of their advantages and disadvan- 
tages . 

The information content of the natural radiation com- 
ponent is high, but has been more of a hindrance than an 
advantage to passive sensing. Existing data interpreta- 
tion techniques are unable to distinguish among the quan- 
tity of radiation contributions and phenomena which dictate 
the characteristics of the passive radiation component. 
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The active component, however, is dependent upon signifi- 
cantly fewer target characteristics than the passive com- 
ponent, and hence less information can be conveyed. The 
smaller quantity of information obtained by active sensing 
may be interpreted more easily, as the complex nature of 
natural radiation is suppressed as noise. The fewer vari- 
ables and higher signal-to-noise ratio are responsible for 
a simpler data interpretation problem, although the number 
of applications is notably restricted by the smaller in- 
formation content. 

Passive sensing has not been utilized to its full 
capability as a result of the limited ability to interpret 
the resulting data. The lack of a suitable model of the 
total extra-sensory radiation situation is the primary 
limitation on interpretation ability. Combined with knowl- 
edge of sensor capabilities and limitations, such a radia- 
tion model would provide the means for more accurate and 
expanded applications for passive sensing. 

Infrared 

The electromagnetic energy spectrum consists of a 
continuum of frequencies from low frequency radio waves 
(100 Hz) to gamma rays (10 2 0 Hz). It is commonly divided 
into bands which exhibit different characteristics, 
although all radiation obeys similar laws of reflection. 
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refraction, and propagation. The term infrared refers to 
that band of the spectrum which ranges from extremely high 
frequency radio waves (3 * 10 1 1 Hz) to visible frequencies 
(3 x 10 1 4 Hz). In terms of wavelength, the infrared band 
extends from 0.75y to lOOOy. Thus infrared differs from 
other radiation (visible, microwave, etc.) only in wave- 
length and frequency. 

The infrared band of the electromagnetic energy spec- 
trum is subdivided into smaller characteristic regions, 
denoted near IR(NIR, 0.75y - 3y) , middle IR(MIR, 3y - 6y) , 
far I R(FI R, 6y - 15y) , and extreme IR(XIR, 15y - lOOOy) . 
Each subdivision corresponds to a band of frequencies for 
which the Earth's atmosphere is fairly transparent (atmos- 
pheric window) or opaque as in the case of XIR. 

The utility of the infrared portion of the electro- 
magnetic spectrum has increased dramatically since the 
military implications became apparent during World War II. 
A survey of unclassified infrared literature reveals in- 
frared systems applications in industry, medicine, and 
science. Indeed, any detection application is possible if 
a measurable variation in radiation is caused by the tar- 
get property of interest. Hudson [10] has compiled a com- 
plete description of the unclassified state-of-the-art of 
infrared applications. 

A major scientific application occurs in the area of 
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qualitative and quantitative remote sensing of earth re- 
sources. Passive detection and processing of infrared 
radiation has been applied to qualitative detection of 
forest fires, underground coal fires, urban development, 
and sea-ice reconnaissance. Quantitative applications in- 
clude Gulf Stream mapping, ground water studies, air pol- 
lution and turbulence measurements , and thermal mapping of 
the Earth's surface. For objects having temperatures on 
the order of 300 °K (approximate Earth temperature), a 
majority of natural radiation emission occurs in the in- 
frared region of the spectrum; 301 falls in the 8-14 micron 
atmospheric window. Thus passive infrared detection is 
particularly suited to accomplishment of the remote sensing 
miss ion . 

Unclassified infrared literature, however, consists 
primarily of qualitative investigations, and quantitative 
studies in which the data is frequently considered to be 
absolute truth with no attempt to evaluate the actual 
quality of the measurements. Scanner type data is often 
viewed as photography, without consideration of emissivity 
or atmospheric effects. This type of data analysis is of 
marginal utility and often of questionable validity. 

A typical investigation has been reported by Wermund 
[11] in which passive microwave and infrared sensors were 
applied to remote sensing of ground water concentrations 
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in arid terrain. Two different conditions, which were 
considered to be indicative of "a high probability of 
ground water”, were dense vegetation and surface moisture. 
Data interpretation, thus consisted of a qualitative search 
for anomalies in the data. Dark areas on infrared imagery 
were interpreted as cold due to the evaporation of surface 
moisture, while areas of low brightness temperature in the 
microwave data were interpreted as areas of low emissivity, 
as a result of the surface moisture. 

In both cases, the results were oblivious of certain 
aspects of the remote sensing problem. Anomalous areas as 
indicated in the data were not necessarily a result of 
ground water, although for the particular situation (arid 
terrain) other geophysical effects were not likely to be 
responsible for the anomalous behavior. Consideration was 
not given to atmospheric effects, background radiation 
variations, or emissivity variations with no apparent data 
deterioration. The consistent nature of the arid terrain 
is responsible for minimization of these various geophysi- 
cal effects and the results of the study show a high de- 
gree of accuracy in prospecting for ground water. The - 
success may be attributed to the relative absence of de- 
teriorating influences in the arid terrain and to the rela- 
tive ease of detection of the anomalous character of ground 


water indicators. 
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Foster [12] has conducted a study of thermal mixing 
in effluent waters using passive infrared imagery in the 
8- 14y range. The results showed an amazing accuracy in 
quantitative remote measurement of water temperature, 
although no attempt was made to account for any geophysi- 
cal variables. Sensor characteristics were also neglected. 
Implementation of three data reduction techniques resulted 
in average accuracies of one to five percent. 

This study has proven the capability to make gross 
temperature measurements under controlled conditions. The 
bulk of the infrared imagery pictured the effluent canals 
of a power generating station. The cross-sectional area 
and flow rate of each canal was essentially constant and 
surface characteristics were invariant except in the area 
of two cataracts. The simplified nature of the controlled 
experiment is responsible for the accuracy obtained. De- 
tailed analysis of sensory and extra- sens ory variables 
would be required for accurate and general applications. 

Scope of Report 

In the following chapters, the total quality of remote 
sensing measurements is considered. The passive infrared 
sensing system is of primary interest; however, the results 
have a degree of impact on active systems, in addition to 
passive sensors of other regions of the electromagnetic 
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spectrum. 

The state-of-the-art of infrared systems and their per- 
formance is introduced in Chapter II. The infrared system 
is introduced schematically and the various components are 
evaluated in terms of their performance figures. The total 
instrument is then considered as a single unit with an as- 
sociated performance figure. 

Chapter III begins the study of extra-sensoring vari- 
ables and their relation to total system performance. The 
topic of interest is the reflective scattering of electro- 
magnetic radiation from target surfaces. General equations 
are formulated which are extremely complex and difficult to 
use . 

In Chapter IV, generation and emission of thermal ra- 
diation is considered. The nature of internal radiation is 
formulated and the scattering theory of Chapter III is par- 
alleled to describe the transmissive scattering of internal 
radiation to the surroundings. The resulting expressions 
are again very complex. 

Conclusions about the theory and its limitations are 
included in the final chapter (Chapter V). In addition, 
recommendations are made for further investigation into the 
extension and simplification of the results which have been 
shown to conform to known results for the special case of 


a smooth surface. 



14 


CHAPTER II 

INFRARED SYSTEM PERFORMANCE 

The performance of an infrared sensing system is 
necessarily dependent upon the characteristics of the com- 
ponents which comprise the total system. Few publications 
in the open literature have addressed the total instrument 
performance problem, due to security classification of the 
subject matter. In this chapter, the performance of each 
component of the system will be considered individually, 
followed by a discussion of the total instrument system 
performance. A schematic representation of the basic 
elements of the modern infrared sensing system: optics, 

detector, cryogenic cooling (for certain types of detec- 
tors), electronic signal processor, and display is shown 
in Figure II- 1 and provides a convenient arrangement for 
the sections that are to follow. 

Optics 

The purpose of infrared optical systems is to form 
images of real objects, which may be emitting their own 
radiation, reflecting radiation from other sources, or 
both. Any distortion in the image or loss or gain of radi- 
ant power introduced by the optical system represents im- 
perfect performance and should be indicated in a perform- 
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ance figure. 

Factors Affecting Image Clarity 

Two processes contribute to distortion of an optical 
image: diffraction, due to the wave nature of radiant 
energy, and aberrations, caused by the optical materials. 
Aberrations may be controlled by the optical designer, but 
diffraction establishes the ultimate limit in optical per- 
formance. Diffraction limited optical systems are cur- 
rently available. 

Diffraction . Diffraction in optical systems occurs 
at the edges of optical elements and at the aperture which 
limits the collected radiation. Even in the absence of 
aberrations, diffraction still causes a point source to 
be imaged as a bright central disk surrounded by several 
alternately light and dark rings. Eighty-four percent of 
the radiant power is contained in the central disk called 
the Airy disk. The angular diameter, 6, of the Airy disk 
is a function of the wavelength of the incoming radiation 
and of the diameter of the collecting aperture, and may 
serve as a partial indicator of the performance of an 
optical system. The relationship is 


8 = 


0.244 X 
' D 0 


(n-i) 
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where: 6 = angular diameter of the Airy disk in millira- 

dians 

A = wavelength in microns 

D 0 = diameter of the collecting aperture in centi- 
meters 

Aberrations . In optical systems, seven types of aber- 
rations can exist, which may be categorized as either 
chromatic or monochromatic. Chromatic aberrations, caused 
by spectral variation of the indices of refraction of the 
lens materials, are of two types: (1) longitudinal, in- 

dicating a variation in the position of the focal point as 
a function of wavelength, and (2) lateral, indicating a 
variation in the size of the image as a function of wave- 
length. Although careful design of refractive optics can 
minimize chromatic aberrations, they are completely elimi- 
nated by the use of reflective optics. 

Monochromatic aberrations, however, may occur in both 
reflective and refractive optical systems and include five 
distinct types: (1) Spherical Aberration, in which rays 

passing through the optics at various distances from the 
optical axis are not brought to a common focus; (2) Coma, 
in which objects located away from the optical axis image 
as an enlarged assymetrical comet shaped blur; (3) Astigma- 
tism, indicated by a point located away from the optical 
axis being imaged as a pair of orthogonal lines; (4) Curva- 
ture of Field, in which the image of a plane object lies on 
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a curved rather than on a plane surface; (5) Distortion, 
in which straight lines not passing through the center of 
the field are imaged as curves. 

Until recently, typical infrared optical systems have 
employed reflective optics predominantly, since relatively 
few optical materials were transparent to infrared fre- 
quencies. Lenses made of ordinary glass which are used in 
the visible range are useful only in the very near in- 
frared, which has provided the impetus for research into 
the infrared transmission of other materials. At present, 
a suitable selection of infrared transmitting materials is 
available to provide a sufficient variety of refractive 
indices and dispersions (variations of refractive index 
with wavelength) so that achromatized compound lenses can 
be constructed. However, all such lens systems have only 
a fairly narrow operating wavelength range. For wide-band 
operation, extending into the middle and far infrared, 
lens systems must be replaced with reflective optics, 
which have focusing properties completely independent of 
wavelength, or perfectly achromatic. 

Although all concave mirrors produce distorted images 
when the target rays arrive from directions deviating con- 
siderably from the optical axis, infrared optical systems 
generally cover only a small instantaneous field of view. 
Also, the detector is normally located near the optical 
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axis. As a result, only certain aberrations are of impor- 
tance to most infrared systems, including spherical aber- 
ration, coma, and chromatic aberrations (if lenses are 
employed) . Aspheric reflective surfaces can completely 
eliminate spherical aberration, although several thin 
lenses also can be combined to reduce it below the diffrac- 
tion limit. The chromatic aberrations may be eliminated 
by the use of reflective optics so that coma is the only 
distortion evident for small field angles and its effect 
decreases with the magnitude of the field angle. 

Factors Affecting Image Brightness 

Radiative power lost or gained in optical processing 
reduces the performance of infrared systems. Diffraction 
can cause a part of the collected radiative power to be 
redirected and subsequently lost. A measure of the dif- 
fraction effect is provided by the diameter of the diffrac- 
tion limited Airy disk, since 841 of incoming radiative 
power is contained within the disk. The linear diameter 
of the disk is related to the angular diameter (II-l) by 
the equivalent focal length of the optics 


(n-2) 


d = Sf 
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where: d = linear diameter of the Airy disk in microns 

f = equivalent focal length of the optics in cen- 
timeters 

6 = angular diameter of the Airy disk in millira- 
dians given by (II-l) 

If the detector is large enough to receive the entire disk 
but small enough to exclude the first diffraction ring, 

16% of the radiation is lost due to diffraction. 

The reflection, refraction, and absorption properties 
of the optical materials can also affect the performance 
of an optical system. For refractive optical elements, 
reflection and absorption represent energy lost; for re- 
flective elements, transmission claims a finite amount of 
radiative power. The spectral transmission of the most 
useful infrared refractive optical materials is shown in 
Figures (II- 2). In general, the transmittance shows a 
distinct spectral variation and does not exceed 90% at any 
wavelength for a 2 mm sample. For the far infrared region 
the transmittance does not exceed 70% for any of the 
materials; at least 30% is lost by reflection and absorp- 
tion . 

Reflective optical elements are inherently more effi- 
cient than refractive elements but performance is also in- 
hibited to some degree by imperfect reflection. The re- 
flectance of aluminum, the most common reflective optical 
material, exceeds 90% for wavelengths longer than 0.5 
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micron and 98% for wavelengths longer than 8 microns, 
another reason for the widespread use of reflective optics 
in infrared systems. 

The reflection from refractive optical elements may 
be theoretically eliminated for a single wavelength by 
covering the surface of the element with an antiref lective 
coating of a specific thickness and refractive index. 
Although the exact thickness and refractive index are im- 
possible to obtain, a 10% deviation in the optimum refrac- 
tive index results in a decrease of only 1% in the effi- 
ciency of the coating. The reduction over a band of wave- 
lengths, such as an atmospheric window, can be significant. 

The emission of radiation by the optical elements 
also deteriorates the quality of the measurement since 
extra non-target raidation is introduced in this manner. 

The emittance of common optical materials is highly clas- 
sified but a reasonable estimate (within 20%) may be 
obtained from Kirchhoff's Law, equating emittance to 
absorptance . 

Optical modulators, called reticles or choppers, are 
sometimes included to provide directional information for 
tracking or to suppress unwanted signals from backgrounds. 
Reticles consist of a pattern of alternately clear and 
opaque areas mounted on a suitably transparent substrate 
and can increase the detectability of a target in the 
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presence of extraneous background detail. The chopping 
frequency is of significant importance to the overall per- 
formance of the system as is indicated in subsequent sec- 
tions of this chapter. 

Modern reflective optical systems for the infrared 
have evolved from the classical reflective systems de- 
veloped for astronomy. Although in a majority of these 
systems the aperture is partially blocked by a secondary 
mirror, the angular resolution is apparently increased 
since a central obstruction reduces the size of the Airy 
disk. Unfortunately the consequence' is a decreased flux 
in the central disk. Popular systems for infrared appli- 
cations include the classical Cassegrainian, Dall-Kirkam, 
and catadioptric systems which employ thin correcting 
lenses (Schmidt, Maksutov, and Mangin Mirror) . A discus- 
sion of modern optical systems for the infrared may be 
found in [14 ] . 

For infrared radiometers and scanners, the fundamental 
purpose is to measure a quantity of radiation incident 
from the target. Therefore the primary optical performance 
factor concerns the percentage of incident target radiation 
which arrives at the detector and the percentage of radia- 
tion arriving at the detector which did not originate at 
the target (factors affecting image brightness) . Using 
the various spectral reflectances, transmittances , and 
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emittances of the optical elements, it is possible to 
arrive at an accurate estimate of the radiative power lost 
and gained. 


Detectors 

The word "detector" has a variety of connotations, 
but may be best described as a transducer of radiant 
energy. In accordance with this definition, there are 
precisely two types of detectors: (1) imaging, which 

simultaneously, form a total image of the target (examples 
are: photographic film or the human eye), and (2) ele- 

mental, which respond to the average irradiance from a 
small field of view (an effective point on the target) . 

It is possible to build an image of a target with elemen- 
tal detectors by sequentially scanning the scene or assem- 
bling a number of elemental detectors into an array thus 
forming an imaging detector. 

It is often more convenient to group detectors into 
two classes that distinguish among the physical mechanisms 
involved in the detection process. Thermal detectors 
operate on the principle that the heating effect of inci- 
dent radiation alters the electrical properties of the 
detector. Response is proportional to the energy absorbed. 
The interaction of incident photons with the electrons in 
a detector forms the basis for the operation of photon or 
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quantum detectors , for which, the response is proportional 
to the number of photons absorbed. Consequently, a plot 
of the theoretical response of photon detectors as a func- 
tion of wavelength for a constant flux is a linearly in- 
creasing function, while thermal detector response is con- 
stant (Figure II-3). 

Performance Criteria 

Past research concerning infrared system quality has 
been largely confined to evaluation of the quality of the 
detector. Several performance figures have been intro- 
duced which describe the partial performance of radiation 
detectors, but a single number representing the total 

i 

performing ability does not exist. 

The detector output per unit input is of primary im- 
portance in calculations of the overall gain function of 
the system. The responsivity of a detector which is 

analagous to its gain is given by 

i 

» = hV 0-31 


where: V s = r.m.s. output voltage of the detector at the 

fundamental chopping frequency 

H = irradiance on the detector 

= detector area 



Photon 



Wavelength 


Figure II-3 

Detector Spectral Response 
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The response time of a detector is represented by its 
responsive time constant which is defined as the time re- 
quired for the detector output to reach 631 of its final 
value after an instantaneous change in irradiance. Many 
detectors respond exponentially, in which case the time 
constant and responsivity are related by 




= ( I + Air's tT 


(n-4) 


where: (R f = responsivity at the chopping frequency f 

Ro = responsivity at zero (very low) chopping fre- 
quency 

x = responsive time constant of the detector 

Although responsivity is useful for gain considera- 
tions, it gives no indication of the minimum detectable 
radiant flux. A measure of minimum detecting ability is 
provided by the noise equivalent power (NEP) . It is de- 
fined as the incident radiant power necessary to produce 
an output signal for which the signal - to-noise ratio (SNR) 
is unity, and is given by 


m c p - HA(jV n _ V n 
V s " R 


(n-5) 
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where: V = r.m.s. value of the noise voltage at the out- 

n put of the detector 

Any actual measurement of NEP must include the electrical 
bandwidth of the noise measurement circuit since the quan- 
tity of noise measured is dependent upon that parameter. 

A similar quantity is the noise equivalent irradiance 
(NEI) , given by 


- H V n _ NEP 
V s A d 

which expresses the irradiance required for a SNR of unity. 
The NEI concept is sometimes used as a measure of the total 
system quality . 

Frequently, for convenience, the reciprocal of NEP, 
known as detectivity (D) , is used as a performance figure. 
The detectivity provides a measure of performance which 
increases (rather than decreases) with the quality of the 
detector, a desirable characteristic in many performance 
factors . 

Detectivity is a function of several variable param- 
eters: incident radiation wavelength, detector tempera- 

ture, chopping frequency, detector bias current, detector 
area, and bandwidth of the circuit used to measure the 




detector noise. 


There is no apparent theoretical relation- 
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ship between detectivity and wavelength or temperature, so 
the variation with these parameters must be found experi- 
mentally. It is customary to make spectral measurements 
at any of several easily achievable and readily reproduci- 
ble temperatures such as those of approximate earth tem- 
perature (300°K), solid carbon dioxide (195°K), liquid 
nitrogen (77°K), and liquid helium (4. 2°K) . The effect of 
chopping frequency is negligible at low frequencies, 

(f << 1/2ttt) and in many cases varies exponentially as in- 
dicated in (II-4). The detectivity has a single valued 
maximum with respect to bias current and is normally ex- 
pressed for that value of current. Extensive theoretical 
and experimental studies have shown that it is reasonable 
to assume that detectivity varies inversely as the square 
root of the detector area, or 

lA 

D A d = constant (IF 7) 


provided the length of the detector sensitive area does 
not exceed five times the width. On the assumption that 
noise voltage per cycle of bandwidth is independent of 
frequency, detectivity varies inversely with the square 
root of the electrical bandwidth Af. 

An attempt to eliminate two of these variables in a 



performance figure has been made by Jones [15] who intro- 
duced the quantity D*: the detectivity referred to an 

electrical bandwidth of 1 Hz and a detector area of 1 cm 2 . 
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D* = D(A d Af)' /2 (H-8) 


This expression is valid under the assumptions of the in- 
verse variation of detectivity with the square root of the 
electrical bandwidth, Af, and the detector area, A^. 
"Dee-star" as it is called, provides a convenient compari- 
son for detectors with different areas used in circuits 
having different electrical bandwidths . The dimensions of 
D* are occasionally referred to as a "jones". The dee- 
star for a specific detector is expressed as D* followed 
in parenthesis by the temperature, chopping frequency, 
and/or wavelength at which the measurements were made. A 
comparison of D* for many common detectors is shown in 
Figure II-4. 

If the performance of a detector is ultimately limited 
by background noise, it is said to be background limited. 
Under these conditions D* is a function of the angular 
field of view of the optical system, and another detec- 
tivity may be introduced which is referred to as an angular 
field of view of tt steradians. Dee- double-star, D**, is 
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Comparison of D* for Common Detectors [16] 
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defined 

D ** = (#f 2 D* tH-9) 

where: ft = effective field of view 

Thermal Detectors 

Infrared thermal detectors, most of which do not 
require cooling, operate by virtue of the heating effect 
of incident radiation on the detector characteristics. 
Theoretically, thermal detectors respond equally to all 
wavelengths, making them ideal for use in radiometers, 
although typical responsive time constants range from a 
few milliseconds to several seconds. 

The first useful infrared thermal detector, the 
thermocouple, was not developed until about thirty years 
following Herschel's experiment. The voltage developed 
at its terminals is proportional to its increase in tem- 
perature. Typically, the electrical resistance is low 
(1 to 10ft); thus, it is well suited for use with transistor 
amplifiers. It is, however, very fragile, generally 
limiting it to use in the laboratory, and its relatively 
long time constant restricts the allowable chopping fre- 
quency to less than 10 Hz. An approximate relationship 
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for the time constant has been formulated by Sanderson 
and stated in [17], 



(n-io) 


where: C = thermal capacity of the thermoelectric junction 

and receiver assembly 

A = energy loss rate of the assembly 
Similarly the responsivity is approximately 



(n-ii) 


Several thermocouples may be connected in series to 
form a thermopile, which displays increased responsivity, 
since the voltages are additive. Increased detector 
resistance simplifies impedance matching to an amplifier, 
but the responsive time constant, usually approaching 
several seconds in length, eliminates the possibility of 
using a chopper. 

A detector which changes electrical resistance when 
heated by radiation is called a bolometer. Thermistor 
bolometers, generally rugged and uncooled, provide good 
field service, displaying high electrical resistance and 
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relatively short time constant (1 to 50 msec). Other types 
of bolometers include superconducting (operating in the 
superconducting regions of certain metals and semiconduc- 
tors), carbon (with D* approximately 10 times that of a 
thermistor) , and germanium (with D* approximately 100 times 
that of a thermistor). 

An optimistic engineering estimate of the minimum 
obtainable NEP with thermocouples or bolometers is provided 
by "Haven's limit". The minimum detectable energy by a 
detector with area of 1 cm 2 is 

E = 3* I O'" joule 


which corresponds to 


NEP 


_ 3 x 10 


HI 


r 


and the maximum D* obtainable is 


Dmax = 1.67 x to 10 T* (H-12) 
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The pneumatic or Golay detector operates on the prin- 
ciple that the pressure of an ideal gas increases propor- 
tionally with temperature increase due to heating by in- 
cident radiation. 

The constructional materials for thermal detectors are 
not in themselves good absorbers and hence must be black- 
ened by applying an absorbant coating. Ideally the coating 
has a uniformly high absorptance at all wavelengths, negli- 
gible thermal capacity, high thermal conductivity, and no 
adverse effect on the electrical properties of the detec- 
tor . 

Quantum Detectors 

The interaction of photons with the electrons in cer- 
tain materials, known as the "photo effect", forms the 
basis for a certain class of detectors. These photon or 
quantum detectors, differ from thermal detectors in that 
they respond to the number of incident photons rather than 
the quantity of incident energy, and hence their spectral 
response is inherently wavelength dependent. Most quantum 
detectors must be cooled to cryogenic temperatures. Their 
performance is generally very good with detectivities one 
or two orders of magnitude greater than typical values for 
thermal detectors, and very short responsive time con- 
stants (usually only a few microseconds). 
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Intrinsic photon detectors are long wavelength limited 
to less than 7 \i since all currently available intrinsic 
materials have forbidden energy gaps in excess of 0.18 
at room temperature, although a slight increase in cutoff 
wavelength can be obtained by cooling. It is also possible 
to control the value of the cutoff wavelength by doping an 
alloy of silicon and germanium. A photoconducti ve detector 
has been produced in this manner and is useful in the 8-14 
micron range, although a minimum of cooling is required. 

The "photo effect" can be classified into two types; 
external, typified by electrons escaping the surface of the 
material, and internal, involving the generation of charge 
carriers within the material. The external type is rarely 
used in the infrared since it has a long wavelength limit 
dictated by the work function of the surface. Photon de- 
tectors based on internal photo effect include photocon- 
ductive, photovoltaic, and photoelectromagnetic type 
detectors . 

Photo con du.cti ve detectors generate electron-hole 
pairs when illuminated by incident photons, resulting in 
a variation in conductivity. The varying conductivity 
modulates the current in a load resistor which is normally 
connected in series with the detector and a dc power 
supply. Any noise in the power supply will deteriorate 
the resulting signal. The responsivity of an extrinsic 
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photoconductive detector of 1 cm 2 area is given by 



v tc/xe 

X cr 


(IT-13) 


where: n 


y 

e 

x 

a 


quantum efficiency (number of charge carriers 
produced per incident photon) 

average carrier lifetime 

mobility of the carriers 

electronic charge 

detector thickness 

detector conductivity 


The maximum theoretical D* for the blip (background limited 
photodetector) condition is 


D*= 2.52 * to' 8 



(n- 14) 


where is the background flux. This theoretical limit is 
included in Figure II-4, p. 31, to show the degree to which 
modern detectors attain theoretical perfection. 

To insure proper utilization of the entire volume of 
the detector, fairly transparent materials are used with 
reflective coatings on the side opposite the illumination 
direction. Intrinsic materials such as silicon, germanium, 
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lead sulfide, lead selenide, indium arsenide, and indium 
antimonide are well suited for photoconductive detection 
as well as germanium doped with gold, cadmium, mercury, 
copper, or zinc. 

A self generating detector (not requiring bias) may 
be constructed by forming a p-n junction in an appropriate 
semiconductor material. Incident photons generate elec- 
tron-hole pairs which are subsequently separated by the 
electric field at the junction, increasing the voltage 
across the junction, hence the term photovoltaic detector. 
The theoretical maximum detectivity is 40 percent greater 
than for phtoconductive detectors, and is given by 


D* = 3.56 x I0 18 



(n- 1 5 ) 


which is also shown in Figure II-4, p.31, although response 
is again limited to 7 y. Common photovoltaic detectors are 
constructed of silicon, indium arsenide, or indium anti- 
monide and are occasionally referred to as "photodiodes". 

Electron-hole pairs generated by incident photons in 
an intrinsic semiconductor material may be separated by an 
externally applied magnetic field. Detectors of this type 
are called photoelectromagnetic and provide very short 
responsive time constants, although the detectivity is 
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lower than that of typical photoconductive or photovoltaic 
detectors . Response to 7 y can be obtained without cooling 
If the performance of a detector is ultimately limited 
by its own internal noise, an increased D* can be obtained 
by placing it in optical contact with a lens. The SNR and 
D* are increased by a factor n g , the refractive index of 
the lens . 

Noise 


Spurious fluctuations (noise) in the detector output 
signal, a result of a stochastic or random process, impose 
limitations on the performance of infrared systems. 

Several types of noise are of importance in infrared 
systems, including Johnson noise, generation-recombination 
noise, radiation or photon noise, temperature noise, and 
several types which vary inversely with frequency, termed 
1/f noise. Two of these, Johnson and 1/f noise, occur in 
both thermal and quantum detectors, while temperature 
noise occurs only in thermal detectors and generation- 
recombination and photon noise occur only in quantum 
detectors . 

Johnson, or thermal noise, occurs in all conducting 
materials: a consequence of the random electron motion. 

The mean square noise voltage, or current per unit band- 
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width, is independent of frequency, a characteristic 
prompting the term "white" noise, and is given by 


V 2 = 4kTRAf (H- 1 6) 


where: k = Boltzmann's constant 

T = detector temperature 
R = detector resistance 
Af = electrical bandwidth of the circuit 

Johnson noise is the predominant noise source in detectors 
at very high frequencies. At slightly lower frequencies, 
1/f noise is apparent. The power spectrum in this region 
is generally expressed, 

p (f) = C (Id ; ] (It- 17) 

where: c = proportionality constant 

I^ c = dc bias current 

The symbols, a and 8 are characteristics of the particular 
detector. Generation-recombination noise is restricted 
to photon detectors and is a result of the random nature 
of e lectron- hole generation and recombination. It is the 
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major source of noise at intermediate frequencies below 
which 1/f noise again prevails. 

The bandwidth of a noise signal is difficult to define 
if it does not display "white" characteristics. To 
alleviate this problem, an equivalent noise bandwidth is 
defined by 


Af = ' gfe - (n ~ l8) 

where: g(f Q ) = P ea ^ value of the noise power gain 

g(f) = power gain at frequency f (Figure II-5) 

Statistically, noise is described by a normal 
(Gaussian) probability density function 

P(x)=^^r exp^jl (H-19) 

where: x = mean noise level (voltage or current) 

a = standard deviation from the mean 

The noise contributed by a single system component 
is represented by its noise factor, F, which is the ratio 
of the SNR at its input to the SNR at its output. It is 
sometimes convenient to express the noise factor in 
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Figure 1 1 - 5 

I'.quivalent Noise Bandwidth 
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decibels, in which case the correct term is "noise figure". 

No satisfactory data have appeared on the range over 
which the response of a detector is linearly proportional 
to irradiance, although most photon detectors respond 
linearly over a range extending to 4 orders of magnitude 
above NEI. Thermal detectors are superior in this respect. 

Cooling 

The response of certain quantum detectors can be 
extended well beyond the 7 micron uncooled limit with the 
aid of cryogenic refrigeration techniques. Several types 
and classes of cryogenic refrigerators are in common usage, 
although a variety of efficiencies are encountered. A 
useful indicator of the efficiency of a refrigeration 
system is provided by its coefficient of performance (COP), 
which is defined as the ratio of the cooling power of a 
system to the input power from all sources. 

Open cycle refrigerators do not provide any reuse of 
the available coolant. Typical systems include liquid- 
transfer , including a coolant reservoir and feed line to 
the detector, Joule-Thomson , employing a Joule - Thomson 
cryostat, solid-re frigerant 3 placed in good thermal contact 
with the detector, and radiative -trans fer , for use in 
space vehicles. 
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Closed cycle refrigeration offers a more efficient 
usage of the available coolant although involving a 
greater power consumption. The open cycle Joule-Thomson 
system may be combined with a compressor to repressurize 
the gas leaving the cryostat, thus forming an efficient 
closed cycle system (COP = 1/100 for refrigeration temper- 
ature = 77°K) . Other closed cycle systems include Claude 
refrigerators (COP = 1/50 @ 40°K, COP = 1/100 @ 20°K, 

COP = 1/200 @ 4.2°K), and Stirling refrigerators 
(COP = 1/30 @ 77°K) , although microphonics represent a 
serious noise problem in both types. 

In the published literature, little consideration 
has been given to the effect of the refrigeration system 
on total system performance. Oscillations or variations 
in temperature, which may be inherent in the cooling system 
or caused by variations in irradiance or detector bias 
current, alter the detector characteristics. Also, noise 
is often introduced by the operation of the refrigeration 
device. For a Stirling refrigerator the noise generated 
when using a Mercury doped Germanium detector is at least 
double that generated by the detector alone (for chopping 
frequencies below 400 Hz or above 2000Hz) . 

Signal Processing 


The output from any detector is a low level signal 
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which must be amplified before it becomes suitable for 
display or control of a recording device. A small pre- 
amplifier is located near the detector to provide a 
suitably high signal level for transmission along shielded 
cable to the remainder of the signal processing elec- 
tronics. Additional amplification, bandwidth limiting, 
and detection are necessary before display is possible. 

For multielement detectors, time multiplexing circuitry 
is also required. Conventional AGC circuitry, logarithmic 
amplifiers, or AGC connected to detector bias circuitry to 
reduce respons ivity , are used to provide the necessary 
wide dynamic range required by many systems to accomodate 

4 6 

changes in signal level on the order of 10 or 10 . 

The most critical element of the signal processing 
equipment is the preamplifier. It must have a low noise 
figure and a minimum voltage gain of 10 to insure that the 
noise contributed by succeeding stages is negligible; a 
gain of 30 to 100 is desirable. AGC is not included in 
the preamplifier in order to maintain the lowest possible 
noise figure. An output impedance on the order of 100 to 
1000Q minimizes pickup in low-noise shielded cables. 

Preamplifiers using vacuum tubes introduce Johnson 
and shot noise to the detector signal. Although maximum 
power transfer occurs when the equivalent resistance of 
the detector circuit (R s ) equals the input resistance of 
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the preamplifier, (R ), a lower noise factor can be obtained 

& 

by prudent selection of R . Assuming 1/f noise is 
negligible, the noise factor is given by 




+ 


Rgq 



(n-20) 


where : 


R 


eq 


R 


R 


equivalent noise resistance of the 
p re amp 1 i f i e r 

source resistance 

grid resistance 


For minimum noise factor, the optimum value for the source 
resistance is 


R s 


opt ^9 


R 


£1 


V? 


Rj + Req 

and the corresponding noise factor is 


(n-2i) 


F • = 

1 min 


4- 2Re< i + 2 

Rn 


Rea 


L Rg 

l Rg/J 


~\/p 


(n- 22) 



47 


The equivalent noise resistance of a triode above 1000 Hz 
is 



2.5 

9 m 


(11-23) 


where: g m = grid-plate transconductance 

For pentodes , 




+ 


81 c 

9m p 


tH-24) 


where: Ic = screen grid current 

p,t = subscripts indicating corresponding 

quantities for the tube connected as a 
pentode or a triode, respectively 


Little information is available in the unclassified 
literature concerning the performance of transistors in 
infrared applications. The noise factor increases at low 
frequencies due to 1/f noise but is constant above 500 Hz 
for high gain planar transistors. In general, a low noise 
factor also requires a small base spreading resistance, 
short circuit common base current gain a approximately 
equal to unity, high short circuit current gain 3 dB 
cutoff frequency, and low emitter current. The optimum 
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source resistance ranges from 100ft to 10,000ft and results 
in a noise figure near 1 dB over a wide frequency range. 

System Performance 

Exact expressions for the performance of composite 
infrared systems are, in general, very difficult to evalu- 
ate but sufficiently accurate engineering approximations 
have been formulated. The most popular approximations 
involve the replacement of wavelength dependent quantities 
with their average over the spectral region of interest. 

For example, the generalized range equation for an infrared 
system limited by detector noise, and viewing an unresolved 
target (that is, one that does not fill the instantaneous 
field of view) has been given by Hudson [18] 



~TrDo(NA)D*JtT X ~ 

2i«Af r (y Vn ) J 



(n-25) 


where: D q = diameter of the collecting aperture 

(NA) = numerical aperture of the optics 

J = radiant intensity of the target over 
the spectral bandpass of the sensor 

t = average transmittance of the atmosphere 
over the spectral bandpass 

x = average transmittance of the optics over 
the spectral bandpass 
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An idealized maximum 
unity for which case 
the noise equivalent 


range, R q , occurs 
the irradiance at 
irradiance NEI, gi 


when V s /V n equals 
the detector becomes 
ven by 


NEI 


2(o 

jAf) 

Vz 

-ttD 0 (NA) 

tfX 


(II- 26) 


Similarly, the idealized maximum range for a background 
limited infrared system employing a photoconductive 
detector and a cooled shield to limit the detector view 
to a half angle 0 is 



~7T Dp J% x S X 

. 4hc 



(n- 27 ) 


where: S = relative detector spectral response 

averaged over the system spectral bandpass 

Dq = detector quantum efficiency 
= background flux 

The generalized range equation can in principle be 
applied to any type of infrared system although certain 
modifications may be necessary. In this thesis, emphasis 
is on passive infrared systems with applications in remote 
sensing of environment: radiometers and line scan imagers. 

For a radiometer viewing an unresolved target, the 
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apparent temperature is related to the system parameters 
by 



2R 2 (coAf)^ 

Lo-A s D 0 (NA) D*7a X 




(n- 28 ) 


where : 


V = peak value of the signal voltage 

A = area of the source 
s 


Although this expression does not give a number value for 
the performance of a radiometer, it does indicate the 
effect of various system parameters. Similarly, the 
apparent temperature of an extended source is given by 



La-Do(NA) D 



14 


(n-29) 


Line scan thermal mapping systems view targets which 
fill the instantaneous field of view at all times and the 
irradiance at the collecting aperature is independent of 
the distance to the target. System performance can be 
expressed in terms of a noise equivalent differential 
temperature (NEAT) which is defined as the temperature 
differential between two adjacent target elements that 
produces a signal equal to the system noise. (The 
temperature differential is actually an apparent 
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temperature differential.) The noise equivalent irradiance 
given by 


NEI = 


2 

D 0 (NA)D* X 


ttCv] 


(n-30) 


where: = velocity to height ratio for the airborne 

platform 

v = pulse visibility factor defined by 
Genoud [19] 


V = 



T 6 a1 


U3-3I) 


In this expression, V ss is the idealized (no loss) output 
from the signal processor and is the detector dwell time 
expressed by 



m-32) 


where: ft = scan rate 

C = number of detector elements 




0 ) 


instantaneous field of view 
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The noise equivalent differential temperature is 
related to NEI by 


m - 331 


where: AN = differential radiance of the target 

which can be expressed using suitable 
engineering approximations (to, t q , 
and s averaged over spectral bandpass): 


aN=%T oS (-|^)dX (n-34) 


Hudson [20] has compiled a table of values for the 
integral in 11-34 for several atmospheric windows and 
target temperature of 300°K. 

Klein [21] has introduced several performance 
indicators which prove useful in the evaluation of 
passive infrared line scanning systems. He has expanded 
on the specific detectivity D* as defined by Jones to 
refer to the detectivity of a detector viewing a hemis- 
pherical blackbody surround at 295°K. 



in- 3 5) 
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where: g = characteristic of the type of detector 

(g = 2 for photoconducti ve , g = /? for 
photovoltaic) 

A = cutoff wavelength for the detector 
c 

D = detector quantum efficiency 
g 

<j> =photon flux density at the entrance 
° aperture 


Similarly, a specific noise-equivalent irradiance is 
defined 


(NElf 


2[co Af] 

mHNAlT^oT 


m-36) 


where: n e = v p / v s 

Incorporating Genoud’s pulse visibility factor and 
a scanning efficiency q , reflecting coverage redundancies 
and waste motion. 


(NED* = 



7TD 0 (NA)T^*J 


C^vj 


(H-37) 


This exhibits the dependence of system performance on 
scanning rate. 
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Noise equivalent irradiance is related to this 
specific NEI by 

NEI=(NEI)*(|jj {H- 38) 

or 

NFT = — K b K r h [tt_-2q) 

LmT^JLc^vJ 111 3y) 

where: = Q b /Q 0 > and K r is the radiation shielding 

factor 

This expression is valid under three assumptions, 
presuming photon limited noise characteristics: (1) no 

solar reflection, (2) background approximate a blackbody 
at and (3) a cold shield, and a cold filter are employed. 

Klein further introduced a figure of merit defined 
as 



(n— 40 ) 



where MDTD is the minimum detectable temperature 
difference as imposed by system parameters. The figure 
of merit is expressed by 
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=[vwmvc£*i] 





in- 4i) 


where: e = emissivity of the target 

(3W./3T) = temp, rate of change of Planck's 
A Law (A- 1) . 

This expression according to Klein, "clearly identifies 
design parameters and target characteristics that 
determine the total system performance of contemporary 
radiation-shielded thermal mappers". 
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CHAPTER III 

ELECTROMAGNETIC REFLECTION 
FROM TARGET SURFACES 

The characteristics of remote sensing data are indica- 
tive of the total extra-sensory radiation situation, in- 
cluding contributions from three general types of sources: 
target, background, and intervening medium. The contribu- 
tion by each source becomes involved in certain electromag- 
netic phenomena including generation, emission, absorption, 
and reflection. An understanding of this extra-sensory 
situation requires a knowledge of the past phenomenological 
history of the specific components involved. Accordingly, 
a theoretical investigation of the reflective scattering of 
radiation from rough surfaces is presented in this chapter. 

General Scattering Problem 

The scattering of electromagnetic waves from "rough" 
surfaces has been a topic of considerable interest in 
recent years. Significant contributions toward a general 
theory of scattering have been made by many scientists 
and mathematicians, although existing theories fail to 
completely explain experimental results. The complexity 
of existing theories and their inability to account for 
experimental data form the impetus for this chapter. The 
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approach presented is not the traditional method used to 
examine rough surface scattering, but is simplified by 
visualization of the physical situation. The scattering 
mechanism is examined by a method which may be strictly 
classified as a Kirchhoff approach, but is not concerned 
with the process of solving th’e boundary value problem. 
The Fresnel reflection coefficients, which represent the 
plane surface solution of the boundary value problem are 
assumed to be valid over small areas of a rough surface, 
allowing the calculation of directional reflectances in 
terms of the surface characteristics . In Chapter IV, the 
approach is modified to examine the emission from a rough 
surface . 

The general problem associated with rough surface 
scattering is indicated in Figure III-l. A quantity of 
electromagnetic energy, incident from the direction r ^ , 
onto an elemental area of surface dS, is redirected into 
various directions depending on the surface roughness 
characteristics. Traditionally, the ratio of the scat- 


tered electric field into a particular direction r , to the 

s 

inci dent electric field from the direction r\ , is referred 
to as the scattering coefficient, Y s ^> which is expressed 


by 
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Figure I I I - 1 

General Scattering Geometry 
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S 

H_ 
i 

*“V 

where the subscripts V, H indicate vertical or horizontal 
polarization in relation to the mean coordinate system. 

For the scattering coefficients, the first subscript in- 
dicates the polarization of the scattered electric field, 
and the second indicates the polarization of the incident 
electric field; this notation follows that of Peake and 
Barrick [22]. The superscripts s, i, indicate scattered 
and incident fields respectively. In this chapter, ex- 
pressions are derived for these scattering coefficients to 
insure that the present development is compatable with 
other work in the field. 

Rough Surface Scattering Model 



The general scattering problem as outlined, may be 
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examined in terms of the traditional concept of scattering 
coefficients. Each incident or reflected direction is in- 
dicated by its azimuth and zenith angles with respect to 
the mean coordinate system (Figure 1 1 1 - 1 , p. 58). Ad- 
ditionally, parameters delineating the surface charac- 
teristics are required to describe the total redirection 
of incident radiation energy. 

Localized Parameters 

Local surface normal . The element of surface dS is 
presumed to contain an "effective area" oriented such that 
the radiation incident from the direction (0^, 4> ^ ) is 
specularly reflected into the direction (0 S , s ) • The 
effective area oriented in the proper direction is a 
function of the surface roughness and of the incident and 
scattering directions . 

Each incident and scattered directional pair (r^, r g ) 
implies the orientation of the effective surface area and 
its local normal, n. The orientation of the local effec- 
tive normal is described by its zenith and azimuth angles 
(0 R , 4> n ) with respect to the mean coordinate system 
(Figure III-2). The statistical properties of the surface 
may be described in terms of these angles, and other pa- 
rameters of interest are related to the zenith and azimuth 
angles of r^ , r , and n. 
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L ocal incidence angle . The local effective normal 
lies in the plane formed by r^, and r , and bisects the 
angle between them in accordance with Snell’s Law of 
reflection (Figure 1 1 1 - 2 , p.61). The angle, 0 ^ , between 
the local effective normal and either r. or r g is of signi- 
ficance in computation of the Fresnel reflection coef- 
ficients and is one half the angle between r. and r g . Thus, 
determination of the local incidence angle may be accom- 
plished by finding the angle between r^ and r g and dividing 
by two . 

A simple geometrical construction (Figure 1 1 1 - 3) 
provides the basis for determination of the local incidence 
angle. According to the Law of Cosines for oblique spher- 
ical triangles [23], the angle between two radii (r^ and r g 
in this case) is directly related to the angles between 
each one and a third radius (F) . The angle between the 
arcs a' and b', defined by (r^, k) and (r s , k) respec- 
tively on the unit sphere, is exactly equal to the angle 
between the projections of rL and r g on the x-y plane 
(<J> s - (j>p . The relationship is given by 

cos 0 is = cos Bs cos#, + sin 9% sin 8 [ cos(<£ s -^) 

(HI- 2) 

The local incidence angle, 0^ n , is given by one half of 
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Local polarization . Incident radiation which is 
polarized vertical or horizontal in relation to the mean 
coordinate system appears to have vertical and horizontal 
components in relation to the local coordinates. The 
vertical and horizontal incident radiation may be trans- 
formed into locally vertical and horizontal components to 
conform to the definition of the Fresnel reflection coef- 
ficients . 

Since the incident fields must be orthogonal to the 
incident direction of propagation, the locally vertical 
polarization direction can be represented by a simple ro- 
tation, ou, of the mean vertical coordinate about the line 
of incidence (Figure 1 1 1 - 4) . The local polarization 
coordinates are oriented such that the local vertical 
coordinate lies in the local plane of incidence formed by 
r^ and r s . The angle of rotation, a^, is obtainable in a 
manner similar to the method used to find the local inci- 
dence angle; using the Law of Cosines for sides of oblique 
spherical triangles, [24] 


cos a-, = 


COS ds - COS g is COS 9\ 
SIN 0 is SIN 0\ 


(m-3) 


For incident radiation consisting of vertical and 
horizontal polarization components, the locally polarized 




66 


components are related by 


• • • 

Elv = Ey COS CLj + Eji| SIN aj 

• • • 

E lh = E|_j COS CLj “ Ey SIN otj 


(m-4) 


Here, the subscripts LV, LH indicate locally vertical and 
horizontal components respectively. The local scattered 
components into the direction r , assuming validity of the 
tangent plane approximation is 


e lv = P v 
e lh = P H e[ h 


(m-5) 


where p.y, are the Fresnel reflection coefficients given 
by [25] 


_ \ cos & - r i l cos 9 m 

v V, COS 0 Z + 77,COS0 in 

(m-6) 

p = tkCOS^n -^iCOS^ 

H %cos6 in + 17 , cos ^ 
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where n 1 , 2 = - 1 v2 = intrinsic impedance o£ the upper, 

£l ’ 2 lower medium 

0. = i 0. = local incidence angle 

in 2 is 

0 2 = local angle of refraction (Figure III-5) 
and from Snell's law [26] 

(m-7) 

where ni , 2 = /y 1 , 2 £ 1,2 = index of refraction of the 

upper, lower medium 

The local polarization of the reflected components 

t • 

j ! 

is related to the polarization coordinates of the reflected 
direction by a simple rotation, a r , of coordinate axes 
about the line of reflection. This transformation is 
analagous to the rotation of incident polarization coor- 
dinates to obtain the local polarization. Therefore, 

<. 

(m- 

The scattered fields in the direction r , are 



cos ft - cos 9 :, cos ft, 

COS a r = ^ 

SIN ft s SIN ft 
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s s s 

Ey — E|^y COS CL|« ** E|^H SIN CLjr. 

(in- 9) 

s s s 

E h - E^jy SIN OLp E lh COS ot-p 

Combining (III- 4) , (III-5), and (III-9), the scattered 
fields may be written: 

Ey — [Ey COS a, + E^SINa,] COS a,. 

• • 

~P H [Eh COS a, - Ey SIN a,] SIN a r 

(m-io) 

E h = [Ey COS a, +• Eh SIN a,] SIN a r 

+ ft [Eh COS a ; - Ey SIN a,] COS a r 


where the Fresnel reflection coefficients are given by 
(III-6) in conjunction with (III-2) and (III- 7) . 

Scattering coefficients . The traditional scattering 
coefficients as defined by (III-l) are obtainable from 
(III -10) ; 
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= /» COSa f COSa r + ^SINa ; SINa r 
7 hh = SINa, SINa r + £ cOSa, COSa r 
7^ v = P, SINa, COS (V ' ft COSa, SINa r 
7y H = /= COSa-jSINa,. - |3 SINa, COSa r 

where p^, are the Fresnel reflection coefficients and 
a., a are the incident and reflected rotations of polar- 
ization coordinates into the local plane of incidence given 
by (III - 3) and (III-8). If the incident and reflected 
directions are interchanged, the theorem of reciprocity is 
not violated by this formulation. These coefficients 
represent the transmission characteristics of plane ele- 
mental surfaces as a function of their orientation with 
respect to the mean coordinate system; they do not account 
for the relative quantities of effective area available 
for specular reflection for each incident and scattered 

directional pair (r., r ). 

^ i ’ s' 
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Surface Model 

Statistical parameters . The characteristics of the 
surface may be completely described by a surface height 
function of position on the mean coordinate plane; 

Z = Z(x,y) (m-12) 

From this approach, other authors have solved the boundary 
value problem at the surface to obtain expressions for the 
scattered fields under certain assumptions. In the re- 
mainder of this chapter the method of physical optics is 
applied to a general surface described by certain prob- 
ability distribution functions. 

The surface distribution of heights above the mean 
plane implies certain other spatial distributions which 
more readily facilitate the present method of solution. 

It will suffice to say that the function (III-12) is re- 
lated to a spatial distribution of local normals, which is 
also sufficient to completely describe the surface charac- 
ters tics . 

Since the local normals are represented by their 
zenith and azimuth angles, measured in relation to the 
mean coordinate system, their statistical properties may be 
described by the statistics of their representative angles 
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9 , <j> n . The spatial distribution of local normals implies 
a probability density function of the random variables 
0 , <j> . Following the notation of Papoulis, [27] the joint 
density for zenith and azimuth angles of the local normal 
is 



(m- 13 ) 


The marginal density of the azimuth random variable 
<j> n is assumed to be uniform and independent of the zenith 
random variable. Thus since <f> may assume only values in 
the region (0 <_ <[> £ 2 it) ; its' marginal density is given by 

1 *^ s ‘£w (0-<^-27r) (rn-i4) 

and the joint density becomes: 


This density function will serve to describe the 
surface characteristics of various random rough surfaces, 
and is related to the specific factors of interest in for- 
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mulation of the scattering problem. 

Effective area . Each element of actual surface area, 
dS , may be associated with a corresponding projection onto 
the mean coordinate plane, dA. A certain fraction, 3, of 
dS is oriented in the proper position to permit radiation 
from the direction r^ to be redirected into the direction 

r . The effective area for reflection from r- into r , 

s 1 s 

within a solid angle dfi, is 

! 

i 

dA e ff =/3 dS (m-16) 

where 3 is the "coefficient of effective area". The ele- 
ment of actual surface dS is a function of the surface 
roughness, therefore 

dS = S dA (m-17) 

where 6 is the "surface roughness factor". 

The coefficient of effective area, 3, and the surface 
roughness factor, 6, are related to the statistical dis- 
tribution of surface normals. The joint density function 
for the zenith and azimuth random variables defines the 
statistics of the surface normals and is given by (III-15). 
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The zenith distribution of surface normals may be of any 
form, depending only on the surface state and will serve 
to describe the surface normal statistics. 

The coefficient of effective area, 3, is related to 
the probability density for the surface normals. The 
probability that the surface normal is oriented in the 
proper direction, within the limits d0 , d<{> is 


1 


2 IT 


a 


d0 n d 


(m-18) 


In spherical coordinates a solid angle bounded by de, 
d<f> is defined 


dn = siN0 d0d<#> (HI- 19) 


Thus, the probability that the surface normal is oriented 
within dft is 



dn 


27TSIN0 


(m-20) 
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The coefficient of effective area is 

U?n) , (ttt 

^ 27T SIN 0 n ® 21) 

This represents the fraction of total surface oriented 
within a unit solid angle about a particular direction. 

The projection of this element of surface onto the mean 
plane is 

ffl(ffn) COS On 

2 ? t sin 9„ 


The total effect of all orientations is obtained by inte- 
gration over the hemisphere 


dA 

dS 


1__ / fplfli) COS On 
S = I 27TSIN On 


dn 


Rearranging and integrating over the azimuth variable 

1 


S = 


/ ir/ z 

f e (0„) COS 0 n 60 n 


(m-22) 



76 


Combining (III- 16 ) and (III-17), 


dAeff - PS dA 


(m-23) 


Substituting (III-21) and (III-22), 



UOn) 

2tt sin 6}, / f # ft cosfl, dO n 

* / o 


dA 

(in- 24) 


which expresses the quantity of effective area available 
for specular reflection from (6^, (J>^) into ( 0 s , <|> ) as a 
function of the surface parameters. 


Phase Interference 


The present formulation has, until now, neglected the 
phase interference due to different path lengths for the 
incident radiation reflected from different elements of 
effective area. In Appendix B, it is shown that there is 
no time average phase interaction of scattered fields due 
to different path lengths if the incident radiation is 
randomly phased. 

Power Scattering Coefficients 

The specific situation involving the scattering of 
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background radiation may be treated assuming that the in- 
cident background radiation is randomly phased. The inci- 
dent goniometric polarized irradiance may be written using 
Poynting's theorem, 


Ixv=e;*h; 
Ixh = Hi 


(m-25) 


where I,„ and I-. u are the vertical and horizontal compo- 
A V AH 

nents of the incident goniometric irradiance and E*, E*, 
H^, Hj*j are the incident electric and magnetic fields from 
the direction (0^, <f>^) . Of this irradiance a certain 
fraction falls on the effective area and is subsequently 
available for specular reflection into the direction 


(e« 


♦.> 


The fraction is expressed: 

=peff _ COS&, 

X »/~P8 COS0; n Av 

^eff = COS 0j - 

^Ah j3S cos 0 m ^Ah 


(m-26) 


and represents the irradiance from (0^, <f>^) which falls on 
the effective area. 

The effective scattered polarized radiance into 
(0 s , cj> s ) from (0^, (j)^) can be written 
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N Xv = Ev*Hv 

(m-27) 

Nxh = Eh x Hh 

where N-^y and are the vertical and horizontal compo- 

s s s s 

nents of the scattered radiance and E V’ E H’ Hy, are the 
effective scattered fields. 

The effective incident and reflected fields are re- 
lated by the scattering coefficients (III-ll), therefore 
the scattered polarized radiance becomes 

N\V " ^V E V +r HV E H> X ( r w H V + 

(m-28) 



Expanding the cross products and realizing 



gives 



(in- 29) 



where the cross products represent the effective irra- 
diance and are given by (III-26). Substituting (III-26) 
into (III-29) , 

_ COS 3 — f y 2 T. + y 2 T ) 

IN Xv pScos&n \ / vv- L Xv 'hv- l Xh] 

„ , x (HI-30) 

M _ _COS_3 lyZ T + -y 2 T ] 

IN XH ps cos 6 m \ 7 HH - l Xh 'vh - L Xvj 

where (66) is expressed by (III-23) and (III-24) and the 

scattering coefficients are given by (III-ll). 

It is instructive to view the reduction of (III-30) 

for the smooth surface case. The density function for 

local normals f n ,(0 , <S> ) becomes a density function for 

0<J> n n 

zenith angle only f Q ^ w ^ich is an impulse at (0 = 0). 

The coefficient of effective area 6 becomes 



fA) 

SIN 0 n 


and the limit as 0 ^ 0 is found, using L’ Hospitals’ rule 
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P = 1 


The surface roughness factor reduces to 


8=1 


and 


cosft = cos0 in 


The scattering coefficients reduce as follows 

^vv = ^ 

^HH = ft 





Therefore Equation (III -30) becomes for a smooth surface 
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N 


Xv 




N 


\H 




(m-3i) 


which is the correct result. 
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CHAPTER IV 

GENERATION AND EMISSION OF 
THERMAL RADIATION 

The radiation component of primary importance in pas- 
sive sensing is the thermal radiation emission from the 
target. It is the predominant carrier of target informa- 
tion, although existing data interpretation techniques are 
unable to extract a majority of the information conveyed. 
This inability to accurately interpret data characteristics 
in terms of the various extra-sensory variables is largely 
a result of the absence of a theoretical model for the ex- 
tra-sensory situation. In response to this data interpre- 
tation problem, this chapter provides the basis for a gen- 
eral theory of the thermal radiation emission. 

Thermal Radiation Generation 

The blackbody radiation laws presented in Appendix A 
are based on certain assumptions about the physical nature 
or radiating material. The assumptions are specifically 
noted in Planck's original derivation [28] and are presen- 
ted in the course of this section. 

Fundamental Concept 


An attempt to describe the various radiation charac- 
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teristics of materials, logically begins by considering the 
actual source of the radiation emerging from the surface. 
The ultimate source is, of course, the atoms and molecules 
which comprise the material, however, a microscopic for- 
mulation of radiation characteristics is not manageable. 
Therefore, following 'Planck , consider a small volume ele- 
ment of the material which encompasses a sufficient quan- 
tity of elementary radiation sources such that homogeneity 
may be assumed (Figure IV-1). The radiative power emitted 
by the volume element is proportional to the volume and 
uniform in all directions, assuming the material is iso- 
tropic and in a state of thermal equilibrium; thermal 
equilibrium in this instance indicates zero temperature 
gradient throughout the material. 

It naturally follows that the radiation propagates 
throughout the medium, interacting with other volume ele- 
ments, and impinging upon the boundary at the surface 
where reflection and transmission occur. Obviously, the 
reflected component remains within the material, while the 
transmitted component escapes the surface (Figure IV-2). 

A basic premise of this chapter is that the transmitted 
component of the internal radiation is precisely the nat- 
ural radiation emission from the material, as it certainly 


must be. 




t -* ®v!\ ^ ^ 

• w I 


V 
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Figure IV-1 


Generation and Emission of Thermal Radiation 
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Internal Radiation Characteristics 

Application of the fundamental concept to describe the 
emission characteristics of an object requires a knowledge 
of the radiation properties within the material. Although 
these are not known specifically, under the given assump- 
tions certain characteristics of internal radiation may be 
noted. For a homogeneous, isotropic medium in a state of 
thermodynamic equilibrium, the radiation flux at any lo- 
cation may be assumed to have perfectly random direction, 
polarization, and phase. From this assumption, the in- 
ternal radiation incident upon the boundary is of the same 
random nature, and the radiative characteristics above the 
surface are then obtainable in terms of the transmission 
characteristics of the boundary. Alternatively, if the 
radiative characteristics above a surface are known, the 
internal radiation situation may be determined. 

Basic radiation theory is based on theoretical objects 
(blackbodies) whose radiation characteristics are such that 
all radiation incident is absorbed (none reflected) and 
reradiated. They are described as perfect absorbers and 
radiators (Appendix A) . The concept of perfect absorption 
and radiation indicates that the transmission of radiation 
energy across the boundary is complete for any direction or 
polarization. Thus, it is a simple matter to infer the in- 
ternal irradiance characteristics of a blackbody. 
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The radiant emittance of a blackbody is described by- 
Planck's Law (A- 1 ,p . 125) , and is related to blackbody ra- 
diance by (A-9, p.129). Since for a blackbody, the trans- 
mission of radiation is complete for any direction, and sub- 
surface spectral goniometric irradiance at any zenith angle 
equals the spectral radiance at the zenith angle; 

Ix 2 - n ‘4 w x 112-11 

Here, the subscript, 2, indicates subsurface. It should 
be noted that this expression also describes the radiation 
situation throughout the blackbody since it represents the 
spectral goniometric irradiance on any plane area internal 
to the blackbody. 

Blackbody internal spectral goniometric volume density 
of radiation energy is related to the subsurface spectral 
goniometric irradiance by the velocity of propagation in 
the blackbody material 

Q L = -Mx=Wlx GE- 2) 

P 

The internal spectral volume density of radiation energy 
is thus. 
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U b x = 4ir I b x (Iff- 3) 


since independence of direction can be assumed without 
error. Similarly, for a physical material, 


m-4) 

and if thermodynamic equilibrium exists in a homogeneous 
isotropic material, independence of direction is assumed. 
Thus 

U x = 4tt { JTT I x 

exactly as given by Planck [29]. 

An interesting question arises: what is the velocity 
of propagation in a blackbody? The answer is not clear 
since the velocity may have any value approaching and in- 
cluding the speed of light. This is a result of the per- 
fect absorption and emission characteristics of a blackbody. 
Complete transmission of energy at the surface of a black- 


(Iff- 5) 
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body requires that its electrical properties assume the 
values characteristic of the medium in which it is placed. 
Thus a blackbody provides a convenient reference and may be 
assumed to allow any velocity of propagation (< c) . 

The internal radiation density of real materials can 
be related to the internal radiation density of a black- 
body at the same temperature. The spectral emissivity, 
e., of the real material is defined in Appendix A as the 
ratio of the internal volume density of spectral radiation 
energy of a material to that of a blackbody at the same 
temperature ; 

U 

e. = — 

X U 

The emissivity is a quality factor expressing a material's 
ability to approximate blackbody radiation characteristics 
internally, and is related to the possible energy transi- 
tions on the atomic scale for the specific material, in 
addition to the spatial density of potential sources of 
radiation. 

Combining (IV- 1), (IV- 3), (IV-5), and (IV-6) , the 
subsurface spectral goniometric irradiance can be related 
to the spectral radiant emittance of a blackbody whose 

permeability and permittivity exactly equal those of the 


b (32- 6) 

X 
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material, since a blackbody may exhibit any velocity of 
propagation ; 

x X2 = IT K 02-7) 

The blackbody spectral radiant emittance, is given by 
(A-l, p. 125), and may be determined experimentally with 
the aid of the development presented in the remaining sec- 
tions of this chapter. 

Peake's Method 

The angular dependence of radiation emission from a 
general surface has been formulated by Peake [30] on the 
basis of thermodynamic considerations. The rationale is 
straightforward; for a material in thermodynamic equili- 
brium with its surroundings, the radiation power emitted 
equals the radiation power absorbed, according to 
Kirchhoff's Law. The total power absorbed or emitted may 
be written 




PE- 81 


where: P a = absorbed power 

P = emitted power 
e r 
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= incident power 
P r = reflected power 


For a homogeneous, isotropic material in thermodynamic 
equilibrium, Kirchhoff's Law may be generalized for spec- 
tral, polarized, radiance and goniometric irradiance. 
Generalizing (IV-8) gives 




Nxh =i 


\H 



( 32 - 9 ) 


where the subscripts V,H represent vertical or horizontal 
polarization components respectively and the superscripts 
e, i, and r indicate emitted, incident, and reflected ra- 
diance or goniometric irradiance respectively. (The ra- 
diometric terminology is defined in Appendix A) . 

The instantaneous power transferred by electromagnetic 
radiation is given by Poyting's theorem 


P = E x H 

where the units are watts per unit area. Extending this 
relationship to goniometric quantities, the instantaneous 
radiance may be written in terms of the spectral reflected 
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polarized gonioraetric fields; 

^Xv" ^Xnv x ^Xav 

^xh " ^xnH x h xuh 


( 31 - 10 ) 


where the subscript fi indicates goniometric fields. 

Using the Fresnel reflection coefficients, (IV-10) may 
be written 

(12-11) 

Similarly, the incident polarized spectral goniometric 
irradiance components are: 

x ' -pi x H ' 

- L Xv~ tl Xav n xav 

lH- 12) 

Ixh- ^XnH x ^xnH 

Combining (IV-9) , (IV-11), and (IV-12), 


N Xv <“ v ^xav x % ^xav 


^Xh ^^xaH xf H^xaH 


N^v = - f? 2 ) 


XV 
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( 12 - 13 ) 

N Xh ' I' - tf ) 4 

These equations describe the angular dependence of polar- 
ized radiation emission from a smooth surface according to 
Peake. The Fresnel reflection coefficients are given by[31] 

cos ft 2 - cos ft 
v \ cos ft + 17, cos ft 

( 12 - 14 ) 

v, cos ft, -17. cos ft, 

p = — 

H t?* COS ft +17, COS ft, 

where: n i , h 2 = intrinsic i mped ances of the 

two media ( /y/e ) 

6i, 02 = angles of incidence, refraction 
(Figure IV- 3) . 

From Snells' Law [32] 



( 12 - 15 ) 
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Figure IV-3 

lixternal Specular Reflection and Refraction Geometry 
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where n^, are the indices of refraction (/ye). The 
angular dependence may be written using (IV-13) . 

Nxv _ 1 . 2. 4 •>?, cos cos 9 ^ 
iL v Kcosfl + ^cosft) 

( 12 - 16 ) 

N XH =1 . -2 _ 4->?,T7 2 cosg, cosg 2 
I^ H H ( Vi cos^ + \ COS 9 , ) 

These equations express the angular dependence of radiation 
emission from a smooth surface as formulated by Peake. An 
alternative approach is now presented. 

Planckian Method 

The radiation emission from an object is considered 
to be that portion of the internal radiation which escapes 
to the surroundings. The internal spectral goniometric 
irradiance upon the surface is assumed to be independent 
of direction and randomly polarized (Lambertian), for a 
homogeneous isotropic material in thermodynamic equilibrium. 
This does not, however, imply that the radiance emitted 
from the upper side of the surface is Lambertian, although 
this is frequently assumed to be the case. The subsurface 
goniometric irradiance, assuming random polarization is: 



96 


^XV2 ^XH2 _ 2 1x2 (12-17) 

where is the total spectral goniometric irradiance 

from medium 2 (subsurface). 

Using Poyntings theorem, the subsurface goniometric 
irradiance may be related to the subsurface incident 
goniometric fields exactly as in (IV- 12). 

^Xv2 = E-xnv x 

(12-18) 

^-XH2 _ x ^xnH 

A fraction of the subsurface irradiance incident on 
the boundary is transmitted across and becomes the radia- 
tion emission. The emitted radiance may also be written 
using Poynting's theorem and the transmission characteris- 
tics of the boundary 
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N X V = T E V2^XttV X "^HV2 ^XftV 

(11-19) 

^XH ’ ^EH2^XflH X HH2^XftH 

where x represents a transmission coefficient; the first 
subscript indicates electric (E) or magnetic (H) , and the 
second indicates vertical (V) or horizontal (H) polariza- 
tion . 

Combining (IV-17), (IV-18), and (IV-19) 

N e =_L t T I 

INI XV 2 'EV2 HV2 - L X2 

P5ZI- 20) 

M e = J_ t T T 

|xl XH 2 / EH2 HH2 ^XZ 

The total emitted radiance is 

n;-n; v *n;„ ®-2h 

and the transmission coefficients are given by 

_ _ 2 y, cos 

EV2 " ^ COS0, + \COS 8z 
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~ _ 2 \ cos 0 2 

nV 2 ” 17, cos 0, + 7? 2 cos 0 2 

( 12 - 22 ) 

y _ 2 T?, COS 0., 

EH2 7?, COS 4 + 7? 2 COS0, 

-r _ 2 % cos 6j 

HH2~ 17, COS 0 2 + 1?;, COS 0, 


From (IV-20) and (IV-22) it is readily seen that the 
angular dependence of natural radiation emission from a 
smooth surface is expressed. It is convenient to define 
a quantity which represents the ratio of emitted radiance 
to subsurface goniometric irradiance; the directional 
transmittance of the boundary is 

0 

2Nxy 417,172 COS 0, COS 02 
^X 2 ( 17 , COS0, + 17 2 COS 02) 2 



r - 2 Nxh 

'XH2 t 


X2 


417,172 COS 0, COS 0 2 

2 

(i7,COS0 2 + 7? 2 COS0,) 


(12-23) 
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It is interesting to note that this expression is 
exactly that as derived by Peake’s method (IV-16), except 
that (IV-23) expresses the ratio of emission to internal 
thermal radiation, while (IV-16) is the ratio of emission 
to externally incident radiation. The difference arises 
from the specific assumptions of thermal equilibrium; 

Peake has assumed total thermal equilibirum with the sur- 
roundings, while the Planckian method involves only the 
assumption of internal thermal equilibrium. Obviously, 
the condition of internal thermal equilibirum cannot exist 
for long if thermal equilibirum with the surroundings is 
not maintained, however, it is a step in the right direc- 
tion toward a general theory of thermal radiation emission. 

For thermal equilibrium with the surroundings, (IV-23) 
reduces to (IV-16) since 

^XVI ~ ^XV2 

IXH, = -^XH2 

The indication is that (IV-23) is a generalization of 
(IV- 16) . 

The emission from a smooth surface is now obtainable. 
Combining (IV-1), (IV- 7), and (IV-23) 
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The radiance of a physical object with a smooth surface is 
thus expressed in terms of the direction, absolute tempera- 
ture of the material, and other physical parameters. This 
formulation also allows the use of the Rayleigh- Je an approx- 
imation when applicable (AT>> 100 ,000y°K) . 


Emission from a Rough Surface 


The transmission of internal radiation to the sur- 
roundings is perturbed by the existance of irregularities 
at the boundary. The general problem associated with 
transmissive rough surface scattering is indicated in 
Figure IV- 4, and may be treated in a manner analagous to 
reflective scattering (Chapter III). A quantity of elec- 
tromagnetic energy, incident from the direction r^, onto 
an elemental area of surface dS is refracted into various 
directions depending on the surface roughness characteris- 
tics. As in Chapter III, it is convenient to introduce 
transmissive scattering coefficients defined as the ratio 
of the scattered electric or magnetic field into the di- 
rection r , to the incident electric or magnetic field 
from the direction r ^ . 
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Figure IV-4 

Transmissive Scattering Geometry 
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where for the field vectors, the subscripts V, H indicate 
vertical or horizontal polarization in relation to the mean 
coordinate system, and the superscripts s, i denote scat- 
tered or incident fields respectively. For the scattering 
coefficients, the first subscript indicates electric (E) , 
or magnetic (H) scattering coefficient; the second indi- 
cates the polarization of the scattered field; the third 
indicates the polarization of the incident field. 


Localized Parameters 


The incident and scattered directions are indicated 
by their azimuth and zenith angles with respect to inter- 
nal (C') and external (C) mean coordinate systems respect- 
ively as in Figure IV-4,(p. 101). The element of surface 
dS is presumed to contain an "effective area" oriented 



such that the radiation incident from (0^, <f>^) is spec- 
ularly refracted into ( 0 s , <J> s ) . The "effective area" is a 
function of the surface roughness and of the incident and 
scattering directions as in Chapter III. 

Local surface normals . Each incident and scattered 

directional pair (r , r ) implies the orientation of the 

i s 

effective area and two local normals: external, n, and 

internal n ' , which are diametrically opposed. The orienta- 
tion of each is described by its zenith and azimuth angles 
with respect to the appropriate coordinate system (inter- 
nal or external) . As a result of the particular defini- 
tion of the coordinate systems, the zenith and azimuth 
angles describing each normal are numerically equivalent 
and are denoted (0 n , 4> n ) ; see Figure IV-5. 

Local incidence angle . The local effective normal 
lies in the plane formed by ? i and r g . The angle between 
r. and n, denoted 0- , is related to the angle between r 

and n, denoted 0 sn> by Snells' Law of refraction, and is of 
significance in computation of the Fresnel transmission 
coefficients. From Snells' Law [33] 


n 


2 


sin 9\ n = n, SIN (9 sn 


( 32 - 26 ) 




Figure IV-5 

Equivalence of Local Normals 
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where ni and riz are the indices of refraction of upper 
and lower media respectively fFig. IV-3, p. 94). The local 
angles of incidence (0^ ) and refraction ( 6 sn ) are related 
to the zenith and azimuth angles of r ^ , r , and n by the 
Law of Cosines for sides of oblique spherical triangles [ 34 ] 
as in Chapter III; 

COS 6 , in =COS0j COS0 n +SIN0|SIN0 n COs(<#> n -<#?) 

(12-27) 

COS 0 S n = C °S 0 S COS 0 n + SI N 0 S SIN 0 n COS (<£ n - ) 


Local polarization . Incident radiation which is 
polarized vertical or horizontal in relation to the mean 
coordinate system appears to have vertical and horizontal 
components in relation to the local coordinates. By the 
same reasoning as in Chapter III, a simple rotation of the 
incident polarization coordinates about the line of inci- 
dence is sufficient to obtain a locally vertical and hori- 
zontal coordinate system. The coordinates are rotated 
until the vertical coordinate lies in the plane formed by 
r^ and r . This angle of rotation, a^, is obtainable using 
the Law of Cosines for sides of oblique spherical tri- 
angles [ 35 ] ; 
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COS <X| 


COS 0 n - COS 0 in COS 
SIN0 jn SIN0j 


( 32 - 28 ) 


For incident radiation consisting of vertical and hori- 
zontal polarization components, the locally polarized 

components are related by 

• • | 

El V = Ey COS aj + E H SIN ex j 

r~ i r- i r- i 

E L h = L h COS aj - LySINotj 

( 32 - 29 ) 

Hlv = Hy COS a i + SIN a, 

• • • 

Hlh = H^cosaj- HySINctj 


Here, the subscripts LV and LH indicate locally vertical 
and horizontal components respectively. The local scat- 
tered components into the direction r g is 
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LV 
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LH' 


T F 1 

/ EH Cl LH 


( 32 - 30 ) 



where the Fresnel transmission coefficients are given by 
(IV-22), and for this geometry 


e, = a 


sn 



( 32 - 31 ) 


The local scattered components are related to mean scat- 
tered components by a rotation a r of polarization coordi- 
nates about the line of propagation. By the same method 
as was obtained, 


COS a r = 


COS Q n - COS B sn COS &$ 
SIN&n SIN 0 S 


(IE- 32) 


The mean scattered fields are 



j— S i— s 

Ey ” Ej^y COS 

E-H “ ^-LH C0S a r 
H* = H* v COS a r 


r— s 

- t LH SIN 
+ E^v SIN a r 
~ H* h SIN a r 


( 12 - 33 ) 


Hh = ^LH COS a r + H^y SIN a r 

Combining (IV-29), (IV-30) and (IV-33) 

Ey = \y [Ey COS a 5 + SIN a,] COS a r 
-7 ^ h [Eh COS a- - Ey SINa.] SIN a r 

E h = 7^ h [Eh cos a. - Ey SIN CLjJ COS a r 

+ 3^y [Ey C0sa { + E^SIN ail SIN a r 

( 12 - 34 ) 

H y = Tj^y [Hy COS a. + Hh SIN aj cos a r 

- ^h[ H H COS a > “ Hy SIN a.J SIN a r 

Hh = ^ H [ H H C0Sa * ' Hy SIN aj COS a r 

+ 7^iv[^V C0S a i + SIN a i] S,N a r 
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where the Fresnel transmission coefficients are given by 
(IV-22) in conjunction with (IV-8) and (IV-31). 

Scattering coefficients . The transmissive scattering 
coefficients as defined by (IV-25) are obtainable from 
(IV- 34) 

T^HH = COS a i COS a r + SIN ctj SIN d r 

^EVV = ”^EV COSa ? COSa r + SIN a. SIN a r 

j ~ y S I N Ot j COS CL^, COS OLj SIN 

T^hv = T^ v COSa.SINa r ~ T^ ( _ j SINa i COSa r 

CE-35) 

3^|j_j|_| ~ COS CL; COS ot-p + y SIN SIN 0^ 

^H VV = ^| V C 0 Sa i C0Sa r + ^hh S,N<X ‘ SINa r 
~^ VH = 7^ v SINa i C0Sa r “7^ H C0Sa.SINa r 
7 ^ hv - 7j^ v COSajSINa r " 7 ^ h SINaj COSa r 
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where x Ey , x EH , x Hy , x HH are the Fresnel reflection coef- 
ficients given by (IV- 22) and (IV-31); o^, a r , are the 
incident and scattered rotations of polarization coordi- 
nates into the local plane of incidence given by (IV-28) 
and (IV-32). Also, if the incident and refracted direc- 
tions are interchanged, the theorem of reciprocity is not 
violated by this formulation. These coefficients repre- 
sent the transmission characteristics of plane elemental 
surfaces as a function of their orientation with respect 
to the mean coordinate system; they do not account for the 
relative quantities of effective area available for spec- 
ular refraction from each incident and scattered direc- 
tional pair (r E , r ). 

Surface Model 

The surface model presented in the previous chapter 
is applicable to the present problem; the effective area 
and actual area per unit mean surface are given in terms of 
the statistical properties of the surface roughness. The 
results are given again here for convenience 



( 12 - 36 ) 
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( 12 - 37 ) 


where: f„(6) = probability density function for the 

zenith angle of the surface normal 

Phase Interference 

As indicated in Appendix B, under the assumption of 
random incident phase, there are no average phase inter- 
ference effects in the scattered field from the elementary 
area of surface. 


Power Scattering Coefficients 


The power scattering coefficients are now obtainable. 
The incident goniometric polarized irradiance may be 
written using Poynting's theorem, [36] 





I 




( 12 - 38 ) 


where: (I.y, I^ H ) = C ve r t ± ° & 1 , horizontal) components 

At! of the incident goniometric irradiance 

E?}, Eu, H, 1 ,, Hu = incident electric magnetic fields 
“ v l! from the direction (0^, c{>^) 
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Of this irradiance, a certain fraction falls on the effec- 
tive area and is subsequently available for specular re- 
fraction into the direction (6 , <f> ) . The fraction was 
given in Chapter III and is 


C 0 S -r- 

Xv = /?S cos 6 m 


T eff . COsfl T 

XH ~/38cos0 in Xv 


( 12 - 39 ) 


representing the irradiance from (0^, <f>^) which falls on 
effective area. The effective scattered polarized radiance 
is 


N 


Xv 







( 12 - 40 ) 


where: = vert i ca l> horizontal components 

' v of scattered radiance 

Ey, E^, Hy, = effective scattered fields 


The effective incident and refracted fields per unit 
solid angle per unit effective area are related by the 
scattering coefficients (IV- 35), therefore the scattered 
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polarized radiance becomes 
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(12- 41) 
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Expanding the cross products and realizing 


e^h^xh^o 


gives 
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(12-42) 
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Xh 


"^EHh"^HHh(Eh X ^H +- ^EHV^HHV Ey x H v 


where the cross products represent the effective irra- 
diance and are given by (IV-39). Substituting (IV- 39) 
into ( IV- 42) , 

M _ COS^i ["-t* y + -y T' T 1 

Xv ~ P 8 cos6L L EVV HVV Xv EVH HVH Xh -I 

(12-43) 

(VI __COS0i J -p y T 1 

■ Xh ps cos0 in L ehh >hhh ^xh ’ehv hhv - l XvJ 
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where g6 is expressed by (IV- 36) and (IV-32) and the 
scattering coefficients are given by (IV-35) . 

For the smooth surface case 


P = S = 1 , cos#i n =cos 9-, 


T = T 

'evv ev 


T =T 

HVV HV 


T =T 

EHH 'EH 


T = T 

'hhh ’hh 


^EVH ^EHV ^HVH %HV ® 

Therefore equation (IV-43) becomes for a smooth surface 


N X V ^EV %v Ixv 

n xh = T eh T hh Ixh 


OZ-44) 


which is the result to be expected. 
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CHAPTER V 

CONCLUSIONS AND RECOMMENDATIONS 

The interpretation of remotely sensed infrared data 
is inhibited by an inability to account for the radiative 
characteristics of targets. Existing performance figures 
for infrared systems viewing extended targets exhibit the 
dependence of system performance on target radiance charac- 
teristics, but current theories fail to completely predict 
the nature of thermal emission and background reflection 
from targets. 

The theory presented in Chapter III provides the basis 
for an analytical model of rough surface reflective .scat- 
tering characteristics. The resulting expressions for the 
polarized reflection may be written as functions of mate- 
rial properties, surface statistics, and incident and 
viewing directional angles only. Two equations result (one 
for each polarization, vertical or horizontal) which des- 
cribe the radiance into any direction due to the reflection 
of the goniometric irradiance on a small area of surface 
from a single direction. The total radiance due to reflec- 
tion of all incident radiation may be obtained through in- 
tegration of each polarized radiance expression over all 
incident directions. The indicated integration appears 
at best to be very difficult, since the argument is an ex- 



116 


tremely complex ratio of trigonometric expressions. 

Analagous expressions describing the thermal radiation 
emission from extended targets are formulated in Chapter 
IV. Again, expressions in terms of material properties, 
surface statistics, and incident and viewing directions 
only, may be obtained although their complexity severely 
limits their utility. Integration over the incident 
directions results in two expressions (one for each po- 
larization) for the thermally emitted radiance. 

The total radiative characteristics of the target are 
expressed by these results. The radiance in any direction 
is the sum of the reflected radiance (Chapter III) and the 
emitted radiance (Chapter IV) , although several assump- 
tions and limitations must be noted. 

The results indicate a degree of correctness as illus- 
trated in the concluding pages of Chapter III and Chapter 
IV, where the expressions exhibit reduction to known re- 
sults for the special case of a smooth surface. This re- 
duction, in itself, warrants further investigation into 
the possible simplification of the resulting formulation. 

In the formulation of the thermal radiation generation 
it was assumed the material was homogeneous, isotropic, and 
in a state of thermal equilibrium, assumptions which are 
not always valid for the general case. Further study is 
required to determine the nature of internal radiation in 
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the absence of these assumptions. The directional trans- 
mission characteristics of the surface are unaffected by 
these assumptions and hence are valid (within limitations) 
for use with any formulation of internal goniometric irra- 
diance on the surface. 

The limitations on the validity of the reflective and 
transmissive scattering theory arise from failure to con- 
sider certain phenomena which are known to occur during 
scattering from rough surfaces. The local radii of curva- 
ture have a dispersive effect on scattered radiation; also, 
multiple scattering and shadowing effects modify the reflec- 
tion and transmission characteristics of the surface. It 
is believed that the theory can be extended to include 
these effects. 

In remote sensing studies employing infrared measure- 
ments, especially images, the relationship between record- 
ed data and the environmental situation is often assumed to 
be a simple one. The radiative characteristics of the tar- 
get and the performance of the sensor are frequently ig- 
nored. This thesis establishes that, in fact, the relation- 
ship is quite complex involving characteristics of the 
sensor in addition to properties of the target. Hence 
considerably more attention must be given to the effects of 
material composition and surface roughness characteristics 
in order to accurately interpret remote sensing data. 
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APPENDIX A 
RADIATION CONCEPTS 

The purpose of this appendix is to provide a firm 
basis for the theory and terminology presented in the main 
body of this thesis. Radiometric quantities, symbols, and 
units have a tendency to be rather cumbersome, and conse- 
quent improper usage is frequent. Successful communica- 
tion among researchers can be accomplished only if a con- 
centrated effort is made toward the use of precise termi- 
nology . 


Terminology 

The large number of terms available to differentiate 
among radiation as a function of various parameters, 

(angle, surface area, wavelength) is generally responsible 
for the resulting misuse. Also, several different sets of 
terms have been developed to describe the various quanti- 
ties of radiometric interest. To aid proper interpretation 
of this thesis, the pertinent symbols, terms, and defini- 
tions are included in Table A-l. As much as possible, the 
symbols and definitions conform to the recommendations of 
the Working Group on Infrared Backgrounds (WGIRB) [37] and 
are consistent with American Standard Z5 8. 1 . 1- 195 3 . 

Several of these terms and symbols require further 
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discussion. Special attention should be given to the term 
goniometric irradiance, designated by the symbol, I, which 
refers to incident radiation power per unit solid angle per 
unit surface area. This term is not generally included in 
tables of radiometric terminology but is useful in the de- 
velopment presented in Chapter IV. 

Another clarification is required, regarding all radi- 
ometric quantities defined per unit area (radiant emit- 
tance, irradiance, radiance, goniometric irradiance) since 
various definitions of unit area are possible. Frequently, 
for the case of radiance or goniometric irradiance, the 
unit area is defined as a unit projection of the actual 
surface area into the viewing direction. This definition 
is useful as it provides a convenient reference area for 
rough surfaces. Blackbody radiance defined in this manner 
varies as the cosine of the zenith angle, obeying Lambert's 
cosine law, and is, hence, termed Lambertian radiance. The 
unit area may also be defined as a unit area of actual sur- 
face or, as in the case of radiant emittance or irradiance, 
a unit area on the mean surface plane. For the sake of 
simplicity, any reference in this thesis to radiant emit- 
tance or irradiance implies a unit mean surface area and 
any reference to radiance or goniometric irradiance im- 
plies a unit projected area, unless otherwise stated. 

The subscripts X and v indicate quantities per unit 
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wavelength or frequency interval. The superscript, b, in- 
dicates blackbody radiation quantities. The symbol e is 
normally used to denote emissivity, a term which is fre- 
quently misused. Emissivity indicates an intrinsic prop- 
erty of a type of material and must be differentiated from 
emissance which includes the effect of the surface char- 
acteristics of a particular sample of the material. In 
this thesis, the greek letter, e, is used to denote emis- 
sance and e' is used to denote emissivity. 

Blackbody Radiation 

Radiation theory is based on theoretical objects, 
blackbodies, which exhibit perfect emission and absorption 
characteristics. Although naturally occurring materials 
may approach blackbody radiation characteristics, true 
blackbodies do not exist. It is possible, however, to re- 
late the radiation characteristics of natural materials to 
theoretical blackbody radiation laws. 

Planck's Law 

The radiation from a blackbody has been formulated by 
Planck on the basis of quantum energy level considerations, 
and may be expressed as a function of wavelength or fre- 
quency of the emitted radiation. The wavelength expression 
for blackbody spectral radiant emittance is 



125 


Wx = c,x 5 [exp(^)-|] ' (A-l) 

where: Ci = 2-rrhc 2 = 3.7415 x 10 4 watt*y 4 /cm 2 

c 2 = hc/k = 1.43879 x 10 4 y°K 

h = Planck's constant = (6.6256 x 10’ 34 W*sec 2 ) 
c = speed o£ light = (2.99 7925 x 10 1 0 cm/sec) 
k = Boltzmann's constant = (1.38054 x 10’ 2 3 W*sec/°K) 
A = wavelength 
T = absolute temperature 

A graphical representation of blackbody spectral radiant 
emittance (per unit wavelength, A) for various temperatures 
is shown in Figure A-l. The spectral radiant emittance 
(per unit frequency, v) is related to the spectral radiant 
emittance (per unit wavelength) by 

W*dX = (A-2) 

and 

dX = X 2 C iv (A- 3) 


Combining (A-2) and (A- 3) 



Spectral Radiant Emittance (W • cm 



FIGURE A-l. 


Blackbody Radiation Curves 
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W £= WxX 2 c 

or 



Rayleigh- Jean's Law 

An approximate expression for blackbody spectral radi- 
ant emittance is provided by Ray leigh- Jean ' s Law, which 
although formulated empirically before Planck's Law, was 
later found to be a special case. Expanding the exponen- 
tial term of (A- 1) in a power series provides the basis for 
obtaining Ray leigh- Jean ' s Law from Planck's Law. Under 
certain conditions (XT >> 100,000y°K) [38] the second and 
higher order terms of the expansion may be neglected, re- 
sulting in Ray leigh- Jean' s approximation: 

W^-|tx 4 IA-5) 

This expression offers a good approximation at very high 
temperatures or at very long wavelengths but obviously can 
not be correct at short wavelengths, since it predicts in- 
finite power as X approaches zero. 
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Wien's Displacement Law 

A maximum in the blackbody spectral radiant emittance 
curves occurs at progressively shorter wavelengths for in- 
creasing temperature, as indicated by the dashed line in 
Figure A- 1 (p . 126). Differentiating (A-l) with respect 
to wavelength, the maximum spectral radiant emittance (per 
unit wavelength) occurs at 

\b - ^3 

*max ~ ~Y~ (A -6) 

where: C 3 = 2 89 7.8y°K 

The value of the maximum spectral radiant emittance 
is obtained by substituting the wavelength of the maximum, 
(A- 5) , into (A-l) 


= C 4 T 5 

max 

where: C 4 = 1.2862 x 10” 15 watt/cm 2 • y • °K 

Stef an-Bolt zmann ' s Law 



(A -7) 


An expression for the total radiant emittance from a 
blackbody is obtained by integrating (A-l) over all wave- 
lengths or (A-4) over all frequencies. Thus, 
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W b =/WxJX=/w'd^=^T 4 (A- 8) 

o o 

where: o = Ste£an-Bolt zmann constant 

= 5.6697 x 10” 12 watt/cm 2 *°K 4 

Blackbody Radiance 

The radiance from a blackbody is independent of direc- 
tion: a consequence of its Lambertian nature and of the 

particular selection of unit area in the definition. It is 
related to the radiant emittance by 

W b = N b 7T (A- 9) 


(Note that the relationship is not W = N2tt as might be 
reasoned from the fact that there are 2tt steradians in a 
hemisphere.) Equation (A-9) was obtained by integrating 
the radiance over the hemisphere 

W =/ N dfi (A- 10) 

27 r 

while paying special attention to the different definitions 
of unit area for radiance and radiant emittance. 
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Non-b lackbody Radiation 

The radiation from natural objects does not obey 
blackbody radiation laws. It is dependent upon such param- 
eters as temperature, surface roughness, wavelength, or 
material properties. But radiation from natural objects 
cannot exceed blackbody radiation for any temperature and 
wavelength. Most materials that are electrical insulators 
obey Lambert's Law for zenith angles below about 60°, while 
materials that are electrical conductors approximate 
Lambertian radiance within engineering accuracy for zenith 
angles less than about 50°. 

To account for their imperfect nature, natural 
materials are said to have an associated emissivity which 
expresses the degree to which a material approximates a 
blackbody. It is commonly defined as the ratio of the 
material's radiant emittance to that of a blackbody. The 
suffix "ivity", however, indicates an intrinsic property of 
a material, which is not compatable with the previous def- 
inition since that ratio includes surface effects. There- 
fore, for the purpose of this thesis, the term emissance, e, 
is used to indicate the ratio of a material's radiant emit- 
tance to that of a blackbody, 
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and the term emissivity, c ' , is used to indicate the ratio 
of the internal volume density of radiation of a material 
to that of a blackbody. Thus, emissance includes the 
effect of the emissivity of a material in addition to the 
influence of surface characteristics. 

Several problems are still inherent in this approach. 
In general, the emissance is known to vary with wavelength 
and surface roughness. Also, typical measurements are of 
spectral radiance, not radiant emittance, since it is dif- 
ficult to measure total radiation over all frequencies and 
directions. For convenience, a spectral goniometric 
(directional) emissance is defined: 



(A- 1 1) 


Angular and spectral dependencies are not known in 
general, although for certain applications, they may be 
averaged and assumed constant over appropriate angular and 
spectral intervals. Also, emissance is an empirical quan- 
tity, which is difficult to measure accurately. In Chapter 
IV a theoretical approach is presented, which clarifies the 
non-blackbody radiation problem through a new approach. 
Spectral goniometric emissance is theoretically derived in 
terms of other physical quantities, including emissivity, 
thus displaying an increased relationship of radiometric 
data to the physical situation. 
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APPENDIX B 

PHASE INTERFERENCE OF 
SCATTERED RADIATION 

Electromagnetic radiation incident upon rough surfaces 
is scattered into both media by reflection and refraction. 
Assuming that the incident radiation arrives in phase, the 
distribution of surface heights introduces phase modulation 
of the scattered radiation. Consequently, the scattered 
radiation pattern exhibits a definite lobe structure due to 
the phase relationships among the scattered components. 
Examination of this phenomena requires the introduction of 
a joint probability density of surface heights and direc- 
tions of the local normal f_. (0 . d> , z) which describes 

the vertical distribution of each effective area. 

Path Length 

The path difference for radiation reflected from ef- 
fective areas at different heights above the mean surface 
is a straight forward geometrical problem. Figure B-l de- 
picts two incident rays, one reflected from an element of 
effective area at height z' and the other at the mean 
plane. The path difference is 

Q z - Z' (cos G; -*• cos $s) (B-I) 
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The element of effective area at each height is now 
assumed to be uniformly distributed over the mean area 
element, dA. The path difference for radiation reflected 
from different elements of effective area at the same 
height above the mean plane is found as follows. Figure 
B-2 shows two incident rays; one is reflected from an 
effective area at the origin of the mean coordinate system 
and the other is reflected from an effective area located 
a distance, d, away, also on the mean plane. The path 
difference is given by 

Jl d = dl ( C 05 + cos <* 2 ) (b-2) 


where d is the separation of the two elemental areas and 


GCX5 <*, = 5 1 


n 0; co 5 (<£,-&) 


CO .5 <*2 9 -5 i Y\ 0 S C 05 a -4u) 


(e -i) 


where ^ is the angle of orientation of d with respect to 
the x-axis of the mean coordinate plane. The coordinate 
pair (d, <j>^) represent a polar coordinate system in the 
x-y plane. The uniform distribution of effective area 
over dA is accomplished by assigning a uniform distribution 
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to the polar coordinates (d, <J>^) . 

Phase Difference 

The total path length difference due to surface height 
and position variations is given by 

1 = i z + JljL (S-4) 


which corresponds to a phase difference of 

(PS) 


It is apparent for this formulation that if the surface 
height variation is small (<<X) , the phase interaction 
becomes negligible, although the surface roughness 
continues to be a factor in the directional reflectance 
and transmittance characteristics. 

Random Incident Phase 

The particular case of random incident phase is of 
special interest in the examination of thermal emission 
and background reflection. Given a number of periodic 
monochromatic waveforms, each with a specific phase with 
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respect to a reference, random phase does not mean that 
instantaneously there is an equal number of waveforms of 
each phase. It means that the probability of any partic- 
ular phase is precisely equal to the probability of any 
other phase but in the instantaneous case there must be a 
certain distribution of phases (Figure B-3). The instan- 
taneous phase distribution of incident radiation may be 
combined with the distributions for d, z, and <p^ and ( B - 5 ) 
to find the instantaneous phase distribution of scattered 
radiation which in general, has some definable form 
different from the incident instantaneous phase distribu- 
tion. 

t 

Consider a single incident phase component, $ 1 , 
represented by a phasor in Figure B-4; the phase distribu- 
tion after reflection is of some definite form which is 
relatable to the incident phase and the surface charac- 
teristics (Figure B-4a). This relationship cannot favor 
any particular phase, therefore, a second incident phase 

! T * 

component $2 rotated in phase by A$' from $1 (Figure B-4b), 
produces a scattered phase distribution which is given by 
a simple rotation by A$’ of all components of the scat- 

t 

tered phase distribution of $ 1 , (Figure B-4b) . 

Now consider an incident distribution of phases as in 
Figure B-5a; the scattered phase distribution is repre- 
sented in Figure B-5b. A phase rotation A$ ' of the 
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(b) 


Figure B-4 

Relationship of a Single Incident Phase 
Component to the Scattered Phase Distribution 
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incident phase distribution (Figure B-5c) produces a phase 
rotation of the scattered phase distribution (Figure 

B - 5 d) . For a sufficient time length, the incident distri- 
bution of phases is uniform, a result of the combination of 
each instantaneous phase distribution. The resultant of 
the instantaneous scattered phase distribution is conse- 
quently uniform, since for a sufficient time span, the 
probability of any particular phase is the same as the 
probability of any other phase. Thus there is no time 
average phase interaction of scattered fields due to 
different path lengths of radiation if the incident 
radiation is randomly phased. 



